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PREFACE 

THE  aim  of  the  present  volume  is  to  set  forth  a  method  of  treating 
the  geometrical  optics  of  any  given  system  intrinsically,  that  is  to 
say,  without  introducing  any  artificial  scaffolding  and,  therefore, 
without  any  arbitrary  splitting  of  the  entities  involved  in  the 
problem  of  tracing  a  luminous  ray  through  the  several  surfaces 
of  a  system.  The  resultant  formulae,  in  vector  language,  of 
course,  as  the  only  appropriate  means  of  intrinsic  expression,  will 
then  be  free  of  any  unnecessary  elements,  and  therefore,  if  finally 
translated  into  any  set  of  scalar  formulae  for  the  practical  use 
of  the  numerical  computer,  will  contain  no  superfluous,  geometrical 
or  arithmetical,  complications. 

Our  purpose  is  not  to  treat  the  whole  subject  of  geometrical 
optics,  but  exclusively,  or  almost  so,  that  part  of  it  which  is  called 
by  the  short  name  of  "  ray  tracing."  This  is  notoriously  the  most 
laborious  part  of  the  computer's  patient  work,  and  becomes, 
without  question,  a  formidable  task  when  he  has  to  deal  with 
skew  rays  and  non-centred  systems.  The  problem  thus  limited 
can  be  put  shortly  : — Given  the  ray  incident  upon  any  system 
of  lenses,  mirrors,  and  prisms,  find  the  emergent  ray. 

The  advantages  of  the  vectorial  method  of  resolving  it  must 
not  be  judged  by  the  conspicuous  shortness  of  the  resultant  formulae 
alone,  but  also  by  the  simplicity  of  their  deduction,  as  compared 
with  the  usual  method,  and  by  the  facility  of  recalling  the  formulae, 
or  of  reconstructing  them  if  forgotten.  Again,  the  nature  of  the 
proposed  method  is  such  as  to  make  the  help  of  drawings — which, 
especially  in  the  case  of  skew  rays,  become,  in  the  best  standard 
treatises,  exceedingly  complicated — almost  superfluous.  This 
circumstance  will  be  particularly  welcome  to  readers  who  are  not 
endowed  with  strong  visualising  powers.  In  connection  therewith 


vi  PREFACE 

the  often  troublesome  discrimination  of  sign  will  become  almost 
automatic. 

As  already  indicated,  the  mathematical  idiom  to  be  employed 
to  secure  these  advantages  and  to  prevent  the  intrusion  of  foreign 
and  artificial  elements  is  Vector  Algebra.  The,  reader  whose 
acquaintance  with  that  natural  language  of  space-relations  is  but 
slight  need  not,  on  that  account,  be  deterred  from  studying  the 
methods  to  be  described.  For  they  will  involve  only  the  rudiments 
of  vector  algebra  :  namely,  the  addition  and  the  scalar  multiplica- 
tion of  vectors  together  with  an  occasional  application  of  the 
vector  product  of  two  vectors.  A  reader  who  is  entirely  ignorant 
of  the  subject  could,  in  a  few  hours,  acquire  this  amount  of  know- 
ledge from  one  of  the  existing  treatises  on  it.  For  instance,  the 
first  twenty  pages  of  the  author's  Vectorial  Mechanics  (Mac- 
millan,  1913)  contain  all  that  is  required  and  even  more.  It  is 
unnecessary,  perhaps,  to  point  out  that  the  rudiments  of  vector 
addition  and  multiplication,  thus  readily  learnt,  will  be  found 
useful  in  other  connections,  such  as  the  study  of  electro-magnetism 
and  kinematics.  It  is,  however,  permissible  to  urge  that  British 
students,  above  all  others,  should  be  more  widely  acquainted  with 
the  language  of  vectors,  so  powerful  and  yet  comparatively  so  little 
used.  For  it  is  their  own  patrimony.  After  due  acknowledgment 
has  been  paid  to  Grassman  it  remains  true  that  the  vector  method 
originated  in  this  country  through  the  creative  genius  of  the 
great  Hamilton,  and  was  shaped  into  a  practical  instrument  of 
investigation  by  Oliver  Heaviside's  pioneer  work  of  simplification, 
seconded  in  America  by  the  independent  and  equally  successful 
efforts  of  J.  Willard  Gibbs. 

In  order  to  avoid  digressions  in  the  text  of  this  volume  it  will 
be  well  to  give  here  a  few  explanations  of  the  symbols  to  be  used 
and  to  collect  the  very  meagre  number  of  vector  formulae  which 
will  be  required  for  our  purpose.  Vectors  will  be  printed  in 
Clarendon,  small  or  capital  type.  Thus 

A,  r,  etc. 

will  be  vectors.  Their  sizes  or  absolute  values  will  be  denoted 
by  the  corresponding  Italics,  as  A,  r  in  the  above  case.  If  r  =  i, 
then  r  is  called  a  unit  vector.  The  vector  sum  of  two  vectors 
will  be  denoted,  as  usual,  by  A+B.  The  scalar  product  of  two 
vectors  A,  B  will  be  written  simply  AB,  and  their  vector  product, 
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which  will  not  often  be  needed,  VAB.  With  these  symbols  the 
only  formulae  required  will  be  : 

A+B=B+A;     AB=BA;     VAB=-VBA, 
A(B+C)=AB+AC;     VA(B  +0)  =VAB+VAC. 

Similarly  for  the  product  of  two  sums  of  any  number  of  vectors. 
If  y  is  the  angle  contained  between  the  two  vectors  A,  B,  then 

A.E=AB  cosy, 

and  the  size  of  the  vector  product  VAB  is 
|  VAB  |  =AB  sin  y. 

The  scalar  autoproduct  of  a  vector  A,  or  its  "  square,"  will  be 
written  AA=A2  =  ^42. 

If  r  is  a  unit  vector,  it  is  evident  that  r2  =  i.  The  angle  y  between 
two  unit  vectors  a,  b,  is  given  by 

ab  =  cos  y. 

The  mutual  perpendicularity  or  orthogonality  of  any  two 
vectors  A,  B  is  equivalent  to  AB=O,  while  their  parallelism  is 
expressed  by  VAB=O.  Instead  of  the  latter  the  form  A  ||  B  will 
sometimes  be  employed.  The  concurrency  (or  equal  sense)  of 
parallels  and  its  reverse  will,  whenever  needed,  be  expressly 
stated.  Finally,  for  any  A,  B,  C, 

VAVBC  =B(CA)  -  C(AB). 

Scarcely  more  than  what  has  here  been  set  down  will  be  required 
for  following  freely  the  course  of  the  deductions. 

The  main  subject  will,  now  and  then,  be  illustrated  by  examples 
or  exercises.  These  will,  for  the  sake  of  better  discrimination  be- 
tween essentials  and  accessory  matter,  be  printed  in  small  type. 

I  gladly  take  the  opportunity  of  expressing  my  best  thanks 
to  Mr.  Frank  Twyman,  Manager  to  Messrs.  Adam  Hilger,  Ltd., 
for  encouraging  this  little  work,  and  to  my  friend  Prof.  T.  Percy 
Nunn  for  reading  and  revising  the  diction  of  the  MS.  My 
thanks  are  also  due  to  the  Publishers  for  the  care  they  have  bestowed 
on  the  book. 


L.   SILBERSTEIN. 


LONDON,  April  1918. 
RESEARCH  DEPT.,  ADAM  HILGER,  LTD. 
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SIMPLIFIED  METHOD  OF  TRACING  RAYS  THROUGH 
ANY  OPTICAL  SYSTEM. 

1.  The  Fundamental  Reflection  and  Refraction  Laws. 

LET  the  unit  vector  n  represent  the  local  normal  to  the  surface 
separating  two  homogeneous  isotropic  media  whose  refractive 
indices  are  /x  and  //.  Whether  n  is  drawn  towards  the  first  or 
towards  the  second  medium,  is  immaterial,  as  will  be  seen  presently. 
Let  the  direction  of  the  incident  ray  be  given  by  the  unit  vector  r, 
that  of  the  refracted  ray  by  r7,  and  that  of  the  reflected  one  by  r", 
the  latter  two  being  again  unit  vectors.  The  plane  r,  n  is  then 
the  incidence  plane.  The  angles  of  incidence  i,  of  refraction  i', 
and  of  reflection  i",  Fig.  i,  are  given,  according  to  the  very  definition 


FIG.  i. 

/L  </>i  .=  -c 

of  the  scalar  and  the  vector  product,  by 

cos  i  =rn,    cos  i'  =r'n,    cos  i"  =r"n, 

or  by  the  sizes  or  absolute  values  of  the  corresponding  vector 
products,  viz.  sinz=|Vrn|, 

and  similarly  for  the  dashed  angles.     Thus  the  law  of  reflection, 
which  asserts  the  equality  of  i,  i"  and  the  coplanarity  of  r,  n,  r", 
can  be  written,  regard  being  paid  to  the  sense  of  r", 
Vr"n=Vrn,    r//n=rn.* 

*  The  second  equation,  besides  being  obvious,  is  also  a  consequence  of  the  first. 
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Multiply  both  sides  of  the  first  equation  vectorially  by  n.  Then, 
by  the  last  formula  of  the  Preface,  and  remembering  that  n2=i, 
the  vector  form  of  the  law  of  reflection  will  become 

r"=r-2n(nr),  (I) 

giving  the  reflected  ray  i"  in  terms  of  the  incident  ray  r  and  the 
incidence  normal  n. 

Again,  the  refraction  law,  which  asserts   that  //sin  i'=/x  sin  i 
and  that  r,  n,  r'  are  coplanar,  is  tantamount  to  writing 

Vns'=Vns, 

where  s  and  s'  stand  for  /xr  and  //r',  respectively,  and  might, 
therefore,  be  called  directed  light  paths.  The  name  is  of  small 
importance.  Keep  in  mind,  however,  that  the  r's  being  unit 
vectors,  s2=/u.2  and  s'2=/i'2.  Treat  the  last  equation  exactly  as 
in  the  case  of  reflection,  i.e.  multiply  it  vectorially  by  n.  Then  the 
result  will  be  s'=s+n(s'n-sn), 

or,  putting  ^p=s/n  -sn=//  cos  i'  -/*  cos  i 

as  a  convenient  abbreviation,  to  be  used  throughout  the  book, 

s'=s+£n.  (II) 

This  is  the  required  vector  form  of  the  law  of  refraction,*  giving 
the  refracted  ray  s'  in  terms  of  the  incident  s  and  the  normal  n. 
It  is  true  that  the  scalar,  or  ordinary  numerical,  coefficient  g  itself 
contains  the  unknown  vector  s',  but  only  in  the  combination 
s'n  =//  cos  i'.  Now  n,  r  being  given,  i  is  given  ;  whence  i',  and 
therefore  also  g,  can  be  rapidly  calculated  by  means  of  the  relation 
/*'  sin  i'=n  sin  i.  It  is,  therefore,  not  advisable  to  eliminate 
V  formally  from  the  right-hand  member  of  (II),  although  this 
can  easily  be  done.  (Vide  infra,  Section  on  Multiple  Reflection.) 
It  will  be  found  more  convenient  to  leave  the  form  (II)  of  the 
refraction  law  with  the  numerical  coefficient  g  as  it  stands. 

*  Formula  (II)  has  also  been  arrived  at,  in  a  roundabout  way,  by  Sommerfeld 
and  Runge,  Ann.  d.  Physik,  35,  1911,  p.  277,  the  object  of  the  authors  being  to 
give  a  general  differential  basis  to  geometrical  optics.  In  our  case  it  appears 
simply  as  the  immediate  vector  translation  of  the  law  of  refraction  considered  as  a 
given  law.  Cartesian  (three  formulae)  equivalents  of  (II)  are,  of  course,  to  be 
found  also  in  some  of  the  old  works.  Formula  (I)  and  its  consequences  were 
given,  in  dj^dic  form,  in  my  paper  "On  Multiple  Reflexion,"  Phil.  Mag.  for 
November  1916. 
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A  glance  at  (I)  will  suffice  to  see  that  the  law  of  reflection  follows 
from  that  of  refraction,  (II),  by  putting  in  g 

-//,  instead  of  //,  (HI) 

a  well-known  rule.  In  fact,  it  is  enough  for  this  purpose  to  multiply 
(I)  on  both  sides  by  the  refractive  index  /x  of  the  medium  common 
to  r  and  r".  Remembering  the  simple  rule  (III),  we  can  henceforth 
use  (II)  alone  as  the  more  general  formula.  Whenever  a  separating 
surface  is  contemplated  as  a  reflector  write  g  =  -  2sn  =  -  2/x,  cos  i. 
Unless  the  contrary  is  expressly  stated  we  shall  conceive  every 
such  surface  as  a  refracting  one. 

Notice  that  if  in  (II),  and  therefore  also  in  (I),  n  is  replaced  by 
-n,  the  formulae  remain  intact. 

Let  now  the  refracted  ray  impinge  upon  another  refracting 
surface.  Then  the  new  refracted  ray  will  be  obtained  from  the 
ray  s'  as  this  was  from  s,  and  so  on,  for  any  number  of  surfaces. 
Call  nlf  n2,  etc.,  in  general  n*,  the  unit  incidence  normals  of  the 
consecutive  surfaces  (at  the  places  of  incidence-  of  the  ray  in 
question)  and  /*«,  /*/  the  indices  of  the  media  separated  by  the 
*th  surface  (Fig.  2),  so  that  /x^'^^+i. 

I 

etc. 


FIG.  2. 

Further  write,  for  each  surface, 

8*  =  PKTK,     **  =  PK'TK, 

as  above  for  a  single  surface,  and  put  gK=slc'nK-sKnlc.  Then, 
applying  the  law  (II)  to  the  successive  surfaces  and  remembering 
that  s/ =8*4-1  (newly  refracted  =  actually  incident  ray),  and  writing 
sl=s  for  the  original  incident  ray,  the  finally  emergent  ray  will  be 
given  by  s'=s +%KnK,  (IV) 

the  sum  to  be  extended  over  all  surfaces. 

In  the  case  when  all  the  surfaces  are  reflectors  the  last  equation 
assumes  a  remarkably  simple  form,  as  will  be  seen  in  a  later  section, 
and  offers  as  a  final  formula  conspicuous  advantages.  In  the 
case,  however,  of  refracting  or  mixed  systems  (refracto-reflectors) 
we  shall  not  actually  require  the  resultant  (IV),  but  shall  employ 
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the  formula  (II)  successively  at  each  surface,  after  having  developed 
an  equally  convenient  formula  for  the  "  transfer  "  of  the  data 
concerning  the  last  refracted  ray  to  the  next  surface.  This  latter, 
the  transfer  formula  (which  will  replace  and  supersede  *  the  formid- 
able array  of  what  the  German  opticians  call  "  die  Uebergangs- 
formeln  "),  is  an  unavoidable  supplement  of  (II),  the  refraction 
formula.  The  former  will  be  soon  developed  in  one  of  the  fol- 
lowing sections,  again  without  the  "  help  "  (which  is  rather  a 
nuisance  or  obstacle)  of  any  scaffolding  or  framework  of  reference. 

2.  The  Nature  of  the  Data  defining  an  incident  Ray. 

It  will  be  recalled  from  the  rudiments  of  vector  algebra  that  all 
of  our  vectors,  such  as  r  or  s,  are  free  or  non-localized  vectors. 
That  is  to  say,  any  line  segment  of  length  s  and  concurrently 
parallel  to  a  vector  s  drawn  concretely  somewhere,  is,  by  the  very 
definition  of  the  equality  of  two  vectors,  for  all  purposes  of  our 
vector  algebra,  equal  to  s.  In  other  words,  to  give  s  means  :  to  give 
fj.  and  the  direction  of  the  ray  only.  Now,  if  the  refracting  surface 
is  a  plane,  this  is  all  that  is  needed  to  find  the  refracted  ray  a 
(/*',  of  course,  being  given).  For  the  whole  planet  has  but  one 
normal  n.  Thus,  unless  we  desire  to  know,  for  instance,  the  glass 
thickness  traversed  by  a  ray  in  a  prism,  the  localization  of  the 
incident  ray  s  is  immaterial  and  not  needed.  Not  so,  however,  when 
the  refracting  interface  of  the  two  media  is  a  curved,  say  a  spherical, 
surface.  For  then  n  is  a  function  of  position  on  the  surface.  Thus, 
s  giving  the  direction  of  the  incident  ray,  we  require  further  data 
for  its  localization.  Now  such  a  datum  is  obviously  the  place 
on  the  surface  where  it  is  struck  by  the  ray  s.  This  amounts  to  two 
scalar  data  defining  position  on  the  surface,  let  us  say  u,  v,  measured 
along  any  network  of  curvilinear  coordinate  axes  spread  out  over 
the  surface,  as  for  instance,  its  lines  of  curvature.  We  need  by 
no  means  specify  these  u,  v.  Clearly,  "  to  give  the  incident  ray  " 
means  to  give  s  |  and  the  place  (u,  v)  of  its  incidence  upon  the 

*  In  the  sense  that  it  will  be  seen  to  be  applicable  to  cases  to  which  the  usual 
formulae  are  not. 

t  Euclidean  space  being,  of  course,  tacitly  assumed  throughout. 

J  Any  two  scalar  data,  the  setter  of  the  problem  prefers.  We  do  not  count  ^ 
involved  in  B  =  JJ.T.  In  general  we  shall  not  count  data  describing  the  optical 
system,  viz.  refractive  indices,  curvatures,  and  so  on. 
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surface  in  question.  But  this  being  given,  the  normal  n  is  also 
given,  which — being  a  unit  vector — amounts  precisely  to  two 
scalar  data. 

Thus,  whatever  the  nature  of  the  numerical  data  for  the  incident 
ray,  be  they  angles  or  projections  upon  any  lines  or  planes,  or 
distances  from  axes  (artificial  or  intrinsic),  we  are  fully  entitled 
to  say  that  "  incident  ray  given  "  means,  for  a  surface  <rK, 
SK  and  nK  given. 

Such  being  the  case,  formula  (II)  finds  the  refracted  ray  s/. 
Moreover,  it  finds  it  together  with  its  localization  ;  for  the  latter 
ray  starts  from  the  point  of  the  surface  <TK  struck  by  the  incident  ray. 
But  in  order  to  apply  (II)  to  the  next  surface  o\+i  we  must  transform 
the  last  acquired  data  S/=SK+I,  n 
into  SK+I,  n^+i, 

that  is  to  say,  we  must  learn  how  to  replace  n*  by  nK+i.  It  is  this 
which  will  be  given  by  the  transfer  formula.  Manifestly,  having 
once  obtained  SK+I,  nK+i,  we  stand  before  the  surface  crK+1  precisely 
in  the  same  position  and  with  the  same  equipment  as  before  the 
preceding  surface.  Thus  the  process  of  tracing  the  ray  can  be 
continued  indefinitely. 

There  is  no  essential  difficulty  in  developing,  in  correspondingly 
general  terms,  the  transfer  formula  for  a  pair  of  surfaces  of  any 
form.  But  for  all  purposes  actually  met  with  in  technical  optics, 
it  will  be  enough  to  develop  the  formula  for  spherical  surfaces  (of 
any  finite  radii),  and  later  on  (in  a  special  section  on  prismatic 
systems)  also  for  the  sub-case  of  two  planes  or  a  plane  and  a  sphere. 
In  the  more  important  case  of  spheres  proper  the  transfer  formula 
is  indispensable  for  finding  even  only  the  direction  of  the  refracted 
ray,  s'K+i,  while  in  the  case  of  planes  (prisms)  the  ray  can  be  traced 
through  the  Whole  system,  as  far  as  its  direction  is  concerned, 
without  the  "  transfer  formula  "  ;  the  latter  being  required  only 
when  one  is  interested  in  the  thickness  of  the  several  media 
traversed. 

3.  The  Transfer  Formula  for  Spherical  Surfaces. 

It  was  said  in  what  precedes  that,  the  point  (u,  v)  on  a  surface 
being  given,  the  incidence  normal  n  is  thereby  given.  In  general, 
however,  .the  converse  manifestly  does  not  hold  ;  two  or  more  or, 
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as  in  the  case  of  a  plane  surface,  even  oo  2  points  having  parallel 
normals.  But  if  the  refracting  surface  is  a  sphere,  of  centre  0,  then 
we  have  a  one-to-one  correspondence  of  points  (u,  v)  and  normals  n, 
there  being  no  two  points  to  which  the  same  n  belongs.  This 
offers  the  advantageous  possibility  of  considering  the  incident 
ray  to  be  given  directly  by  s  and  n.  It  will  be  convenient  in  the 
case  of  any  spherical  surface  to  draw  all  the  normals  n  as  unit 
vectors  from  the  sphere's  centre  as  their  common  origin.  The 
two  vectors  s,  n  being  given,  we  know  already,  by  (II),  how  to 
trace  the  ray  through  the  corresponding  surface,  i.e.  how  to  find 
the  refracted  ray  s'. 

What  we  still  require  is  to  make  this  ray  ready  for  tracing  it 
through  the  next  surface  of  the  system.  This  is  an  easy  geometrical 
problem,  although  the  optical  authorities  have  made  it  a  compli- 
cated one  by  treating  it  with  too  many  "  auxiliaries."  Treating 
it  (according  to  our  principle)  intrinsically,  it  is,  and  will  remain, 
a  simple  problem. 

Let  o-j,  o-2,  as  representatives  of  any  pair  <rK,  o^+i,  be  two  con- 
secutive refracting  spherical  surfaces  of  the  given  optical  system. 
Let  0lf  02  be  their  centres,  and  Rlt  R2  their  radii.  Draw  all  the 
unit  normals  n^from  Olt  and  all  n2frcm  02,  as  a  rule.  The  optical 
system,  as  always,  being  given,  so  are  Olt  02,  and  all  the  following 
centres,  03,  etc.,  collinear  or  not  (centred  or  not  centred  system). 
Thus  let  D  be  the  scalar  distance  of  the  two  centres  contemplated 
and  a  an  unit  vector  drawn  from  0^  towards  <92,  as  in  Fig.  3.  (The 


FIG.  3. 

reader  should  note  that  the  surfaces  need  not  turn  their  convex 
sides  towards  one  another,  nor  need  0^  precede  02 ;  the  formulae 
to  be  given  presently  will  hold  always.)  Let  nx  be  the  particular 
normal  pointing  to  Plt  the  place  struck  by  the  incident  ray  sx 
which  penetrates  into  the  second  medium  as  BI  =s2,  and  strikes  the 
next  surface  in  P2,  a  point  to  which  the  particular  normal  n2 
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corresponds.     Having  found  S1'=s2  by  (II),  we  wish  to  find  P2, 
i.e.  the  normal  n2. 

In  order  to  obtain  it  walk  round  the  quadrilateral  0-^P-f^O^ 
back  to  the  starting  point  Ov  This  being  a  closed  polygon,  no 
matter  whether  plane  or  not,  or  whether  you  recross  your  path  * 
or  not,  the  vector  sum  of  the  four  vectors  traversed  is  nil,  i.e., 
calling  L  the  (scalar)  length  PiP2,  or  medium  "  thickness  "  traversed 
by  the  ray,  and  remembering  that  S1/=s2=/i2r2,  and  therefore, 


(i) 

This,  in  fact,  is  already  the  required  transfer  formula,  expressing 
the  unknown  n2  in  terms  of  things  given  —  all  but  the  scalar  co- 
efficient L,  a  length.  In  order  to  find  this  (which,  by  the  way,  is 
not  only  needed  for  "  tracing  "  the  ray  but  will  serve  the  computer 
in  his  further  task  of  calculating  "  astigmatism  and  curvature  "), 
eliminate  from  (i)  the  principal  unknown  n2  by  squaring  both  sides 
of  this  vector  equation.  In  doing  so  remember  that  n22  =  i,  as 
are  also  the  squares  of  the  remaining  three  Clarendons.  Then  the 
result  will  be 


a  quadratic  equation  for  L,  with 

P=JDa-#1n1=D-R1,  (3) 

which  is  a  given  vector.     This  finds  L  and  completes  the  required 
transfer  formula  (i)  which  can  be  written,  more  shortly, 

£2n2=Lr2-P.  (i') 

Is  not  this  simple  ?  In  order  to  estimate  the  simplification  of 
the  deduction,  due  to  the  vector  method,  read  the  same  subject, 
for  instance,  in  the  well-known  collective  work  of  the  scientific 
Jena  staff,!  where  it  covers  ten  pages  of  print  (in  which,  to  be 
just,  not  one  but  three  alternative  "  Uebergangs  "-methods  are 
described,  none  inferior  to  the  others  either  in  rigour  or  in  clumsi- 
ness). In  harmony  with  this,  their  "  refraction-formulae  "  are 
equally  complicated,  the  whole  apparatus  consisting  of  something 
like  fifteen  or  sixteen  formulae.  The  complication  in  the  scalar 

*  Which  would  happen  with  the  choice  of  another  exemplifying  figure. 
t  Die  Bilderzeugung  in  optischen  Instrumenten,  edited  by  M.  von  Rohr,  Berlin, 
J.  Springer,  1904,  pp.  52-62.     It  is,  nevertheless,  a  valuable  work. 
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treatment  is  of  course  due  to  the  circumstance  that  skew  rays 
are  aimed  at,  not  cutting  0^0^.  But  as  the  reader  will  have 
noticed,  we  have  nowhere  assumed  that  the  incident  ray  has  any 
such  privileged  orientation.  To  repeat,  the  quadrilateral  01P1P202 
need  not  be  plane  ;  the  resulting  formulae  (i)  or  (i')f  and  (2) 
with  (3)  are  perfectly  general. 

The  solution  of  the  quadratic  equation  (2)  will,  in  numerical 
cases,  cost  but  little  trouble,  and  this  will  be  repaid  by  obtaining 
L  not  only  as  an  auxiliary  for  the  transfer  formula,  but  also  as 
one  of  the  "  elements  "  needed  for  estimating  astigmatic  and 
associated  defects  of  the  optical  system  in  question.  The  equation 
(2)  has,  of  course,  in  general  two  roots, 


L  =  (Pr2)  ±  v/(Pr2)2  +  #22-P2  (4) 

corresponding  to  the  two  points,  P2  and  P2'  (of  Fig.  3),  in  which 
the  ray  cuts  the  second  sphere,  and  in  practice  there  will  be  no 
difficulty  in  discriminating  between  the  root  or  point  (P2  in  the 
case  of  Fig.  3)  to  be  retained  and  that  to  be  rejected.  And  if 
both  roots  coincide  or  are  complex  (conjugate),  there  is  an  end 
to  the  ray-tracing  ;  for  in  these  cases  the  ray  does  not  penetrate 
into  the  next  medium. 

It  would  be  entirely  superfluous  to  enter  at  this  stage  into  the 
details  depending  on  the  configuration  of  Olf  O2,  and  the  following 
centres,  together  with  the  sense  of  the  curvatures  (concave,  convex  ; 
convex,  convex,  and  so  on)  .  The  reader  will  do  well  to  draw  several 
figures  for  the  possible  cases,  especially  for  those  in  which  the  sides 
of  the  quadrilateral  cross  one  another,  —  and  to  inspect  them  care- 
fully. Among  other  things  he  will  then  convince  himself  that,  with 
the  vector  a  drawn  always  from  Ol  to  Oa  (whether  Oj  precedes  0,  or 
follows  it),  formula  (i),  supplemented  by  (2),  (3),  remains  always 
valid. 

Having  obtained  n2  by  means  of  the  above  formulae,  we  can 
proceed  to  trace  the  ray  through  the  surface  having  its  centre 
at  02  ;  similarly  we  shall  pass  from  that  to  the  surface  centred 
at  03,  and  so  forth.  No  complications  will  be  met  with  if  it  should 
happen  that  the  system  is  not  centred,  i.e.  if  0lf  02,  03,  04,  etc., 
are  not  collinear,  or  even,  not  coplanar.  For  the  first  two  surfaces, 
a=a12  will  be  the  "  axial  "  unit  vector  ;  for  the  second  and  third 
another  unit  vector  a23,  drawn  from  02  towards  03,  will  replace 
it,  and  so  on.  These  vectors,  being  not  artificial  but  intrinsic, 
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inherent  in  the  optical  system,  cannot  give  rise  to  any  superfluous 
complications  of  the  calculator's  work.  If  the  centres  happen 
to  be  arrayed  in  a  straight  line,  all  the  better  :  in  that  case  a  single 
unit  vector  a  will  perform  the  function  of  the  axis  throughout 
the  tracing  process.  The  calculations  will  become  correspondingly 
simpler,  again  for  an  intrinsic  reason,  to  wit,  because  the  optical 
system  as  such  is  a  particularly  simple  sub-case  of  the  general 
non-centred  system. 

4.  Summary  of  the  Tracing  Process,  for  Systems  of 
Spherical  Surfaces. 

It  will  be  well  to  summarize  the  whole  procedure  of  tracing  a  ray 
through  any  such  system. 

The  incident  ray  being  given  in  direction  and  position,  the  data 
of  the  problem  are,  in  whatever  form, 

SJL    and    n^ 

From  these,  /zx  cos  i^—^^  and  i±  by  j*/  sin  i±  '  —^  sin  i,  follow 
immediately,  whence  the  scalar  coefficient 

£!=/*/  COS  *j'-)"l  COS  *!.* 

This  being  found  we  have,  by  the  refraction  formula  (II), 


which  is  also  s2,  the  ray  incident  upon  the  second  surface  of  the 
system. 

The  vector  drawn  from  Ol  to  the  next  centre  02  being  D12  =  £>12a12, 
and 

P=D12-JR1n1, 

the  medium  thickness  L  is  calculated  as  the  appropriate  root 
L  =  (Pr2)  ± 


of  the  quadratic  equation  (2),  with  r2  = 
Finally,  by  the  transfer  formula  (i)  or  (i'), 

n2  =  ~(Lr2-P). 
"s 

Thus  the  two  vectors  are  found, 

s2     and    n2, 
ready  to  trace  the  ray,  thus  localized,  through  the  second  and 

*  Which  is  to  be  replaced  by  -  2^  cos  ^  ,  when  the  surface  is  contemplated  as  a 
reflector. 
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to  carry  it  over  to  the  third  surface,  and  so  forth,  up  to  the  last 
surface  of  the  system,  giving  the  finally  emergent  ray  through 
two  vectors,  say, 

s'    and    n', 

the  first  determining  its  direction,  and  the  second  localizing  it 
completely  with  respect  to  the  system. 

The  results,  at  any  intermediate  surface  <TK,  as  well  as  at  the 
last  surface  of  the  system,  are  all  of  the  same  form  as  the  original 
data  (sp  n^. 

Each  of  these  pairs,  whether  of  data  or  results,  amounts  to  the 
knowledge  of  four  scalar  numbers,  two  involved  in  s  (not  counting 
/x)  and  two  in  the  unit  vector  n. 

This  exhausts  the  matter  in  its  general  aspect,  the  procedure 
being  valid  for  any  orientation  and  position  of  the  incident  ray 
and  for  any  system  of  spherical  surfaces,  centred  or  not. 

With  regard  to  its  particular  or  practical  aspect,  the  successful 
application  of  this  compact  set  of  vector  formulae  to  numerical 
cases  will  depend  only  upon  one  negative  condition,  which  is  :  not 
to  spoil  their  intrinsic  simplicity  by  inappropriate  splitting  or 
decomposition  of  the  vectors.  A  large  portion  of  the  remainder 
of  the  present  work  will  be  dedicated  to  an  attempt  to  show  the 
reader  how  to  avoid  this  danger. 

As  a  matter  of  fact,  readers  skilled  in  the  handling  of  vectors 
will,  after  what  has  been  said,  hardly  require  any  further  explana- 
tions. This  little  book,  however,  being  chiefly  intended  for  those 
who  are  not  experts  in  the  use  of  vectors,  it  has  seemed  indispensable 
to  add  some  hints  at  least  which  are  likely  to  help  the  majority 
of  readers  in  applying  the  general  formulae  to  more  or  less  concrete 
cases.  These  will  be  given  in  some  of  the  following  sections. 

5.  Hints  and  Examples. 

Having  hitherto,  in  developing  the  fundamental  and  general 
part  of  the  subject,  followed  a  systematic  course,  we  can  henceforth 
without  danger  of  confusion  depart  therefrom  and  choose  our 
subjects  freely.  Our  aim  being  to  explain  the  general  procedure 
by  illustrations  and  hints,  a  pedantic  observance  of  order  would 
only  hamper  us  unnecessarily. 

Let  us  start  with  the  simple  case  of  a  lens  consisting  of  a  single 
piece  of  glass  of  index  /x  limited  by  two  spherical  surfaces  <rlf  o-2, 
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the  surrounding  medium  being  air  or  empty  space.     To  fix  our 
ideas  let  us  take  a 

Biconcave  Lens.  Let,  as  before,  01}  Rl  be  the  centre  and  the 
radius  of  the  first  surface,  upon  which  the  incident  ray  s  =  r 
impinges,  and  02,  R%  the  centre  and  the  radius  of  o-2,  the  back 
surface  ;  finally  let  nlf  n2  be  the  unit  normals,  drawn  from  the 
centres  towards  the  points  at  which  the  ray  enters  and  emerges, 
as  in  Fig.  4,  which  however  is  scarcely  needed.  By  assumption, 


FIG.  4. 

the  two  unit  vectors  r,  n  are  given,  amounting  each  to  two 
numerical  data.  The  first  practical  question  is  :  How  are  they 
given  ?  Generally  speaking,  this  will  depend  upon  the  decision 
of  the  person  who  sets  the  problem.  Yet,  since  we  are  here  setting 
our  own  problems,  let  us  put  the  data  in  a  reasonable  form,  seeking 
only  to  avoid  artificial  frameworks  of  reference.  Now,  the  lens, 
consisting  but  of  two  spherical  surfaces,  offers  us  but  one  intrinsic 
directed  element,  viz.  its  axis,  which  we  will  represent  by  the  unit 
vector  a  drawn  from  Ol  towards  02, — these  points  themselves 
having  already  been  exploited  by  taking  them  as  the  origins  of 
the  normals.  Besides  a  there  is  nothing  intrinsic  to  build  with 
a  framework.*  Now,  if  the  ray  r  happens  to  lie  in  a  meridian 
plane  (passing  through  the  axis),  then  nothing  more  than  a  is 
needed  for  the  ray's  relevant  definition  ;  say  the  angles  r  and  nx 
make  with  the  axis.  But  let  r  be  a  skew  ray,  not  contained  in  any 
meridian  plane. t  Then  the  lens  itself  offers  us  manifestly  nothing 
to  complete  the  determination  of  the  direction  r.  But  the  ray 

*  A  plane  normal  to  a  may  be  a  welcome  auxiliary,  but  gives,  manifestly,  nothing 
independent  of  a  itself.  Nor  shall  we  look  for  the  rim  of  the  lens  to  help  us. 

t  The  simpler  case  of  a  meridian  or  non-skew  ray  may  be  taken  up  as  a  sub-case, 
and  as  an  easy  but  useful  exercise. 
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itself,  as  localized,  does  offer  us  something,  and  exactly  as  much 
as  is  needed.  In  fact,  it  strikes  the  first  surface  in  a  given  point  Plt 
given  in  one  way  or  another.  Through  that  point  draw  your 
'*  first  meridian  "  plane,  marking  it  with  any  unit  vector  b,  say 
perpendicular  to  a,  and  count,  from  that  plane  or  semiplane,  the 
geographical  longitude  <£,  the  companion  of  the  polar  distance 
or  colatitude  ^1  =  ^a,  nx.  Thus  6,  <f>  will  be  a  polar  system  of 
coordinates  with  centre  Ov  Remember  that  we  have  drawn  our 
first  meridian  plane  a,  0lf  b,  only  through  the  end  point  Pl  of  r, 
but  that  the  ray  is  in  general  not  contained  in  that  plane. 

Now,  this  is  all  we  want.  The  vector  a  being  furnished  by  the 
lens,  and  b  by  the  incident  ray,  we  have  everything  which  is 
required  to  write  down  numerically  (scalarly)  the  data.  In  fact, 
"  r  given  "  means  now  ra,  rb  given,  say  : 

ra  =  cosa,    rb  =  cos/2.  (5) 

For  symmetry,  although  it  is  superfluous,  introduce  a  third  unit 
vector  c  =  Vab  (normal  to  a,  b,  towards  the  reader),  and  write 

re  =  cos  y,     so  that  (ra)2  +  (rb)2  +  (re)2  =  I .  (6) 

Keep  in  mind,  however,  that  the  first  two  are  by  themselves  suffi- 
cient. This  settles  the  question  of  r.  Next,  "DJ  given"  now 
means  0lt  fa  given,  i.e.,  to  be  general, 

%  =  a  cos  Ol  +  (b  cos  fa  +  c  sin  fa)  sin  0lt  (7) 

and  since,  by  the  choice  of  the  first  meridian,  fa  =  o, 

n1=acos^1+bsin^1.  (7  a) 

In  fine,  the  complete  numerical  data,  immediately  attackable 
by  the  computer,  are 

flu     <£i=°»     ra>    rb> 

the  latter  two  as  in  (5),  or,  if  he  prefers,  with  re  as  in  (6)  instead 
of  rb. 

Such  being  the  scalar  data,  act  according  to  the  compact  pro- 
cedure of  Section  4.  Thus  * 

cos  z'j  =  rjiij  =  cos  Ol  (ija)  +  sin  6>1  (ijb) 

-  cos  #! .  cos  a  +  sin  0X  cos  /?,  by  (70),  (5) ; 

whence  sini/s-  sinij     and    ^1=ficosi1/-cosi1, 

*  Remembering  that  in  the  present  case  Bl=Tl=Tt  /*}  =  !,  /*1/=/*a  =  Ai,  ^—\. 
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ready  for  the  first  refracted  ray 

P*l=tl+gl*l=^2-  (8) 

If  the  reader  wishes  to  split  this  at  the  present  stage  he  can  do  so 
at  once,  multiplying  (8)  scalarly  by  a,  b  or  c.  But,  being  interested 
only  in  the  finally  emergent  ray,  let  us  proceed  according  to  the 
general  prescription  of  Section  4.  To  "  transfer  "  the  ray  r2, 
found  by  (8),  to  the  surface  cr2,  write  D  for  the  distance  0-f)^  and 
P  =  Z)a-7?1n1,  as  before;  then  the  scalar  product  required  for  L 
will  be 


D  R 

=  -(cosa  +  glcosOJ--i(cosil  +  gJtby  (8),  (5),  (;«), 

and,  the  square  of  P, 

P*  =  D*  +  R1*-2DRlcos  Olt  by  (7  a). 

This  finds  the  glass  thickness  traversed,  which  in  the  present 
case  is,  manifestly,  the  smaller  of  the  two  roots  of  (2), 

L  =  (Pr2)  - 


which  the  reader  will  properly  contract  by  using  the  last  two 
expressions.  The  "  transfer  "  is  thus  completed,  giving  n2,  and 
therefore  also  the  point  P2,  where  the  ray  strikes  the  back  surface, 

na  =  i(Lr2-P). 

To  know  the  numerical  coordinates  of  this  point  may  be  of 
immediate  interest.  Now,  calling  02,  fa  its  colatitude  and  its 
longitude,  with  02  as  centre,  but  still  using  the  previous  first 
meridian  plane,  in  short  the  same  unit  vector  b,  we  have  just, 
as  in  (7), 

n2  =  -  a  cos  62  +  (b  cos  fa  +  c  sin  fa)  sin  02) 

i.e.  for  the  required  two  angular  coordinates, 

cos  B2  =  -  n2a,     cos  fa  sin  02  =  n2b,     sin  <j>2  sin  02  =  n2c, 

or,  ultimately, 

cos  B2  =  -  an2  ;     tan  fa  =  ^2.  (9) 

Notice  that  fa  being  zero,  fa  need  not  necessarily  be  so;  in 
other  words,  n2  is  in  general  not  coplanar  with  a,  nA.  It  may  be 
well  to  recall  here  that,  for  any  three  vectors, 

AVBC  =  BVC  A  =  CVAB 
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is  the  volume  of  the  parallelepipedon  constructed  on  A,  B,  C  as 
edges,  and  that,  therefore,  the  sufficient  and  necessary  condition 
of  coplanarity  of  n1?  a,  n2  is  the  vanishing  of  the  (scalar)  expression 

R&Vnpz  =  LaVn1r2  -  aVi^P, 
or,  since  P  is  itself  coplanar  with  a,  nx,  of  the  expression 

//aVn1r2=aVn1r1,  by  (8). 

We  thus  see  that  <£2^0,  i.e.  that  P2  is  outside  the  meridian  plane, 
unless  the  incident  ray  rx  happens  to  be  in  a  meridian  plane.  For 
any  skew  incident  ray  we  have  </>2  =£  o. 

The  coordinates  <£2,  #2  of  the  leaving  point  P2  being  thus  found 
by  (9),  the  emergent  ray  is  localized,  and  it  remains  only  to  find 
its  direction  r2'.  This  is,  by  the  refraction  formula, 


the  latter  to  be  found  at  once  from  the  known 
r2n2=-^(L-Pr2), 

as  was  g,  a  moment  ago.  This  completes  the  tracing  of  the  ray. 
For,  if  a',  /?'  be  the  angles  which  the  finally  emergent  ray  makes 
with  a,  b,  we  have  only  to  multiply  (10)  scalarly  by  a,  b,  to  obtain 

cos  a'  =r2'a,     cos  /3'  =r2'b. 
Thus  the  original  incidence  data  having  been  put  in  terms  of 

*ii     <£i  =  °,     ai     A 

the  final  answer  is  obtained  in  similar  language,  viz. 
6,,,      +»      a',      ff. 

This  essentially  exhausts  the  problem  of  a  simple  lens.  If 
instead  of  a  biconcave  we  have  a  biconvex,  or  a  "  vex-cave  "  lens, 
etc.,  the  only  modifications  needed  will  concern  the  sign  of  a 
in  the  expansion  of  n2,  and  the  sign  of  the  square  root  appearing 
in  the  "  glass  thickness  "  L  traversed.  The  exhaustive  discrimina- 
tion of  all  possible  cases,  which  will  among  other  things  depend  also 
upon  the  relation  of  Rlt  R2,  Dt  may  be  left  as  an  easy,  but  useful, 
exercise  for  the  reader.  The  central  glass  thickness  is,  in  the 
above  case,  D  —  (Rl  +  R^),  and  similarly,  mutatis  mutandis,  for  the 
remaining  cases. 

Any  further  information  concerning  the  emergent  ray  may, 
whenever  required,  be  easily  deduced  from  the  obtained  02,  <£2;  a',  /?', 
or  vectorially,  from  n2,  r2'.  Thus,  for  instance,  to  decide  whether 
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the  emergent  ray  meets  at  all  the  axis  of  the  lens  (at  a  finite  or 
infinite  distance),  test  the  triple  product 

£=-aVr2'n2,  (n) 

which  is  a  scalar,  or  ordinary  number.  If  p  =  0t  the  emergent 
ray  or  its  prolongation  meets  the  axis  somewhere  or  is  parallel 
to  it ;  and  if  p  =£  o,  then  r2'  is  a  skew  ray.  By  (10)  and  (8)  we  have 

i),  and  so  on. 


Let  the  reader  complete  the  development  of  pt  and  discuss  the 
circumstances  of  its  vanishing  or  not.  A  useful  exercise. 

As  a  further  exercise,  assume  (abstractly  or  in  a  numerical  example) 
that  p  ±  o,  so  that  the  emergent  ray  is  a  skew  ray,  and  find,  vectorially, 
by  using  r2',  n2,  the  shortest  distance  of  the  emergent  ray  from  the 
lens  axis  ;  also  the  point  of  the  axis  which  lies  at  that  distance 
from  the  ray. 

As  another  exercise  treat  the  trivial  case  of  p=o  and  find,  again 
vectorially,  the  distance,  from  the  back  surface,  of  the  point  of 
intersection  of  the  ray  ra'  and  the  axis. 

Construct  a  handy  formula  for  the  angle  between  the  original 
incident  ray  and  the  finally  emergent  one,  remembering  that  its 
cosine  is  r^a'. 

The  reader  will  also  find  it  useful  to  invent  for  himself  similar 
problems,  resolve  them  vectorially  and  evaluate  the  appropriately 
scalarized  expressions  numerically.  In  this  way  only  will  he  acquire 
the  needed  skill  in  handling  not  only  vectorial  optics  but  vector 
investigations  in  general,  which,  after  all,  will  be  found  also  to  be  a 
comparatively  pleasant  occupation. 

Not  confining  himself  to  lenses,  let  him  also  treat  a  single  or  two 
spherical  mirrors,  for  which  gK  assumes  the  simpler  form  -  2  cos  iK, 
as  explained  before. 

In  actual  numerical  ray  tracings  the  principal  and  the  accessorial 
vector  formulae  will  best  be  written  out  (and  glanced  at  whenever 
necessary)  in  the  margin  of  the  computer's  sheet,  or  on  an  independent 
little  sheet  of  paper  conveniently  placed  before  his  eyes.  Now  and 
then  it  will  be  found  convenient  to  note  a  vector  result  by  writing 
out  the  figures  found  for  its  "  components,"  as  for  example 

P  =  i5a- 24-75^, 

having  previously  decided  upon  the  choice  of  the  length  unit.  It 
will  be  remembered,  of  course,  that  every  unit  vector  is  dimensionless, 
being  a  vector  divided  by  its  own  (scalar)  size  or  representative 
length.  Finally,  details  concerning  the  manner  of  writing  and 
arranging  the  numbers  or  their  logarithms  will  best  be  left  to  the 
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individuality  of  each  reader.  It  is  a  familiar  experience  that,  within 
very  wide  limits,  everybody  finds  his  best  way  of  arranging  the 
numbers  by  himself,  after  he  has  repeatedly  executed  a  given  type 
of  computation.  It  is  by  no  means  desirable  to  drill  away  human 
individuality. 

6.  Prismatic  Systems. 

Leaving  for  the  present  the  spherical  separation  surfaces,  let  us 
turn  to  plane  ones.  A  series  of  different  optical  media  separated 
from  one  another  by  plane  surfaces  a-K  constitutes  a  prism — or 
prismatic  system.  This  is  of  almost  equal  interest  and  importance 
with  lens  systems.  As  a  matter  of  fact,  any  prism  system  is  but 
the  simplest  sub-case  of  what  has  been  treated  in  the  first  two 
sections.  Yet,  in  order  not  to  omit  any  opportunity  of  familiarizing 
the  reader  with  the  vector  treatment  of  optical  problems,  let  us  enter 
somewhat  into  the  details  of  these  comparatively  simple  systems. 

A  single  Prism.  Let  us  begin  with  the  case  of  but  two  plane 
surfaces  <rlt  cr2  limiting  (to  an  extent)  a  medium,  say  glass,  of 
refractive  index  /*.  Outside  /*  =  i.  In  the  first  place,  let  us  treat 
the  rays  apart  from  their  localization,  i.e.  only  as  far  as  their 
directions  are  concerned.  Each  refracting  face  is  then  fully 
characterized  by  a  single  unit  normal,  nt  for  the  first,  n2  for  the 
second  plane  face.  Remembering  that  inversion  of  the  n's  leaves 
all  the  fundamental  formulae  intact,  draw  both  normals  towards 
the  air,  as  in  Fig.  5,  so  that  the  prism  angle  or  refracting  angle  e 
is  given  by 

HjUg  =  COS (TT  -  «)  =  -  COS  e.  (l2) 

The  question  of  an  "  auxiliary  "  framework  now  entirely  dis- 
appears. Unlike  the  simple  lens  (or  any  centred  lens  system)  the 
prism  is  liberal  enough  to  furnish  us,  by  its  own  intrinsic  elements, 
a  complete  natural  reference  system.  For  it  has,  together  with  nlf  n2, 
a  well-determined  edge.  Take  a  unit  vector  e  along  the  edge, 
normally  to  the  plane  of  Fig.  5,  and  towards  the  page  (or  down- 
wards if  your  reading  table  is  horizontal),  so  that 

Vn1n2  =  esin  (TT  - «)  =  e  sin  e.  (13) 

The  a.ngle  c  being  given,  of  course,  this  is  a  simple  relation  between 
the  three  natural  reference  vectors.  Returning  to  Fig.  5,  drop 
the  outline  of  the  prism  itself,  shift  the  normals  n1}  n2  to  a  common 
origin,  and  draw  instead  the  unit  edge  e,  as  in  Fig.  50.  The  prism 
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will  then  be  represented  by  the  simple  three-vector  skeleton,  with 
/x  subintended. 

Now  let  any  ray  r  impinge  upon  (the  "  first  "  face  of)  the  prism 
thus  represented.     This  given  ray  may  also  be  drawn  from  the 


FIG.  5.  FIG.  5*. 

same  origin,  though  it  is  scarcely  necessary.  It  finally  emerges 
from  the  prism  as  r2'=r',  say,  viz.  by  (IV),  Section  I, 

i'=T+gi*i+gp*  (14) 

Here  r  is,  of  course,  in  general  any  skew  ray.  The  explicit 
expressions  for  the  scalar  coefficients  gl}  gz  can  be  found  at  once. 
But  even  without  the  knowledge  of  their  values  the  vector  formula 
(14)  shows  us  at  a  glance  an  interesting  and  perfectly  general 
property  of  the  prism.  Multiply  it  scalarly  by  e  ;  then  the  result 

will  be 

r'e=re=cos0,  say.  (i$) 

Thus  the  final  emergent  ray  makes  with  the  edge  the  same  angle  0 
as  the  original  incident  ray,  whatever  its  orientation.  In  other 
words,  the  only  effect  of  the  prism  is  to  turn  the  ray  rigidly  round 
the  prism  edge  as  axis.  The  amount  of  this  rotation  depends  upon 
the  properties  of  the  prism,  ^  and  e,  and  upon  the  orientation  of 
the  incident  ray. 

The  ray  r  being  given,  so  is  the  incidence  angle  i  =  arccos  (rn^, 
whence  also  the  coefficient  glt  appearing  in  the  resultant  formula 
(14), 


Similarly  the  second  coefficient  gz  could  be  written  down,  as  in  the 
case  of  the  lens  (and  somewhat  easier).  But  it  is  more  elegant  to 
represent  fully  the  final  emergent  ray  r'  without  the  help  of  #2-* 

*  Especially  as  the  orderly  way  of  calculating  g^  would  presuppose  the  know- 
ledge of  the  final  exit  angle  *'  =  z'2'. 
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To  get  rid  of  this  troublesome  coefficient  multiply  (14)  vectorially 
by  n2 ;   then,  by  (13),  and  remembering  that  Vn2n2  =  o  identically, 

(17) 


sin  e 

This  expresses  a  vectorial  property  of  the  finally  emergent  ray 
in  terms  of  things  which  are  all  known.  Together  with  the  capital 
property  (15)  the  last  equation  completely  determines  the  emergent 
ray  in  its  orientation. 

Exercises.  Discuss  equation  (17),  deducing  from  it  various 
vectorial  consequences  and  translating  each  of  these  into  geometrical 
language. 

Find  the  expression  for  the  amount  of  rotation  round  prism  edge 
as  axis,  which  converts  the  incident  r  into  the  emergent  r'. 

Find,  in  the  general  case  of  any  skew  r,  the  limit  condition,  beyond 
which  the  ray  does  not  emerge  but  is  internally  reflected  (at  o-2). 

Treat  vectorially  the  trivial,  but  instructive,  case  of  r  _L  e-  (Assume 
at  the  outset  r  of  the  form  ant  +  bn2  ',  a,  b  scalars.)  Translate  your 
results  into  scalar  language,  deriving  thus  the  usual  formulae. 

In  practice  (spectroscopy,  etc.)  the  aim  is,  of  course,  to  make 
re=o,  or  r_Le.  Yet  the  manufacturer,  and  the  student  in  his 
physical  laboratory,  can  but  approach  this  ideal  within,  say,  ±  10* 
or  even  broader  limits.  It  is,  therefore,  important  to  estimate 
the  effect  of  this  error  upon  the  emergent  r'.  Develop,  therefore, 
the  above  formulae  up  to  the  very  end  by  assuming  a  value  of 

re  =  cos  0 
slightly  differing  from  zero,  i.e.  equal  to  a  small  fraction.     In  other 

words,  put  #=-±8,  where  8  is  a  small  angle,  say  lo*  or  i7  or  5'. 
Employ,  whenever  necessary,  the  well-known  developments  in  series  : 

Xs      X6 

smx=x — J-  +  -J-..., 

^2    5x\  (18) 

~2!  +  4l~""J 

Find  the  effect  upon  r7  of  a  small  change  or  "  variation  "  in  the 
direction  of  the  incident  ray  r-  Replace  r  by  r-f  8r,  remembering 
that  ra  =  i,  and  therefore,  r8r=o,  or  8r_|_r,  and  so  on. 

Similarly  study  the  effect  of  a  slight  variation  8e  in  the  refracting 
angle  of  the  prism,  as  expressed  by  (12).  Scalarize  the  resulting 
formulae  and  interpret  them  in  physical  (geometrical)  language. 

Medium  Thickness  traversed.  In  what  precedes,  the  rays 
incident  upon,  transient  through  and  emerging  from  the  prism 
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were  considered  with  regard  to  their  direction  only,  independently 
of  their  localization.  In  fact,  it  was  possible  for  that  purpose  to 
replace  the  prism  by  a  mere  triad  of  unit  vectors,  n^  n2,  e.  The 
"  transfer  formula  "  was  not  needed  owing  to  the  circumstance 
that  the  same  n2  served  as  incidence  normal  for  the  whole  second 
plane.  This  is  the  reason  why  the  treatment  of  a  prism  is  con- 
spicuously simpler  than  that  of  a  lens.  In  short,  the  tracing  of 
rays  through  a  prism  as  far  as  their  direction  goes  requires  no 
14  transfer  "  process.  Some  device  of  that  kind,  however,  becomes 
necessary  if  we  desire  to  find  the  length  of  the  ray,  whose  direction 
is  r/  =r2,  within  the  prism,  or,  as  we  will  call  it  for  brevity,  the 
medium  thickness  traversed.  Let  this  be  L. 


FIG.  6. 


Obviously,  in  order  to  find  it,  we  must  localize  the  incident  ray  r 
with  respect  to  the  prism,  by  giving,  say,  the  point  Px  where  it 
strikes  the  first  face  o-j.  Since  Px  cannot  be  defined  by  rij  (as  in 
the  case  of  a  spherical  interface),  the  normal  being  common  to  all 
points  of  the  plane,  it  is  indispensable  to  determine  position  on  o^ 
by  some  coordinates.  As  such  we  may  conveniently  take  ylt 
measured  along  e  (the  edge),  from  an  arbitrary  origin  0  on  the 
edge,  and  x^  perpendicularly  to  the  edge.  Similarly  y2,  xz  f°r  the 
second  face  o-2  of  the  prism.  We  can  collect  these  pairs  of  co- 
ordinates by  introducing  the  vectors 

(19) 


where  k1?  kg  are  two  unit  vectors  along  the  axes  of  xlt  x2,  as  in 
Fig.  6.    ' 

To  obtain  L,  walk  round  the  triangle  OP-JP^G  ;    then  the  result 
will  be 

P2=p1+Lr2.  (20) 
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This,  in  the  present  case,  is  the  "  transfer  formula,"  somewhat 
simpler  than  in  the  case  of  spherical  surfaces.  Multiply  (20)  scalarly 
by  n2,  which  is  normal  to  p2 ;  then 

Zjyij^-p^,  (21) 

which  finds  L,  and  completes  at  the  same  time  the  prismatic 
transfer  formula.  Also,  of  course,  L2  =  (px  -  p2)2. 

Exercises.  In  (21)  substitute  TZ=TI,  using  the  fundamental 
refraction  formula,  thus  completing  the  little  investigation. 

Study  L  in  its  dependence  upon  the  prism  angle  e,  and  upon  the 
direction  and  the  position  of  the  original  incident  ray. 

Deduce  the  formulae  for  the  simplest  sub-case  re  =  o. 

Find  the  deviation  of  the  ray,  noting  that  its  cosine  is  rr'.  Establish 
the  conditions  of  minimum  deviation. 

System  of  any  Number  of  Prisms.  The  tracing  of  any  ray, 
in  direction  and  position,  by  the  successive  application  of  the 

refraction  formula  / 

s*+i  =SK  =&K+gKn.K 

and  of  the  transfer  formula  typified  by  (20)  with  (21),  offers  no 
difficulties,  nor  essentially  novel  points,  and  need  not,  therefore, 
detain  us  here.  We  can  imagine  that  by  the  application  of  the  said 
process  all  the  scalar  coefficients  gK  are  already  determined  by  the 
given  incident  ray  r  and  by  the  properties  of  the  given  system 
of  prisms. 

Thus,  supposing,  for  instance,  the  first  and  the  last  medium  to 
be  air,  the  emergent  ray  r'  will  be,  by  (IV), 

r'=r+%n,c=r+£1n1+  ...  +gKnK,  (22) 

where  n*  is  the  unit  normal  of  Kth  plane  face.  In  this  more 
general  case  (*>2)  it  will  be  convenient  to  draw  all  the  normals 
nx,  n2,  n3  and  so  on  in  one  and  the  same  sense,  say,  from  the  first 
towards  the  second,  from  the  second  towards  the  third  medium, 
and  so  on.  Again,  as  in  the  case  of  a  single  prism,  the  whole 
prismatic  system  can  be  represented  by  a  skeleton  of  coinitial 
unit  normals  nlf  n2,  n3,  etc.,  and  unit  edges  e12,  623,  etc.  These 
latter,  no  matter  whether  coinciding  in  space  or  not,  can  also  be 
shifted  to  the  same  origin.  (An  air  gap  will  also  count  as  a  com- 
ponent "prism,"  with  ^  =  i.)  Thus  a  system  of  N-i  prisms, 
i.e.  of  N  interfaces,  will  be  entirely  represented  by  a  skeleton  of 
N  vectors  n*  and  N  -  i  edge  vectors  e12,  etc.  If  elf  €2,  etc.,  are  the 
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prism  angles,   the  relationship  of  the  latter  to  the  former  unit 
vectors  will  be  expressed  by 

sin  ex  .  e12  =  Vn^g,     sin  e2 .  e23  =  Vn2n3,  etc. 

Any  given  incident  ray  r  impinging  upon  this  system  will  be 
converted  into  r',  as  expressed  by  (22),  where  the  gK  are  assumed 
to  have  already  been  found.  (The  system  itself  constitutes,  of 
course,  the  best  and  more  than  complete  framework  of  reference, 
thus  giving  no  place  for  artificial  "  auxiliaries.")  All  properties  of 
the  emergent  ray  can  at  once  be  derived  from  the  vector  equation 
(22).  We  know  already  that  the  sole  effect  of  each  prism  taken 
separately  is  to  turn  the  incident  ray  round  the  edge  as  axis, 
leaving  the  angle  0  with  the  edge  intact.  The  effect  of  the  whole 
system  will  be  the  resultant  of  a  number  of  such  rotations  round 
e12,  e23,  etc.,  as  axes,  taken  in  their  proper  order,  of  course, — if,  as 
in  general  will  be  the  case,  the  edge  vectors  do  not  coincide,  i.e.  if 
the  prism  edges  are  not  all  parallel  to  one  another. 

The  optical  manufacturer  attempts  in  ordinary  circumstances 
(as  in  the  case  of  spectroscopes)  to  make  the  prism  edges  parallel, 
i.e.  e12  =  623  =  634,  etc.  But,  of  course,  he  can  succeed  in  doing  so 
only  approximately.  Let  us  call  such  an  ideal  system  a  correct 
one,  and  any  other,  with  edges  more  or  less  deviating  from 
parallelism,  a  defective  prism  system.*  Now,  especially  if  a  system 
of  prisms  is  defective,  formula  (22)  may  be  very  helpful,  because 
with  HJ,  n2,  n3,  etc.,  not  coplanar,  any  of  the  ordinary  Cartesian 
methods  would  soon  lead  to  the  greatest  confusion  and  would 
require  very  careful  drawings  and  possibly  three-dimensional 
models  to  follow  the  rays  in  their  tortuous  course.  The  vector 
formula  will  therefore  be  particularly  useful  if  it  is  required  to 
estimate  the  effect  of,  say,  small  deviations  from  parallelism  upon 
the  emergent  ray. 

If  the  system  of  prisms  happens  to  be  (or  is  considered  as)  a 
correct  one,  then  all  the  normals  nx  are  coplanar  and  all  the  unit 
edge  vectors  are  fused  into  one,  e,  say.  Multiply  (22)  by  this 
vector  ;  then  the  result  will  be 

r'e=re, 

i.e.  0'  =  0,  precisely  as  for  a  single  prism,  as  could  be  expected. 
Whatever  the  orientation  of  the  incident  ray  r,  the  prism-system 

*  Although  the  reader  may  well  imagine  cases  in  which  non-parallelism  is  a 
feature  deliberately  aimed  at. 
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turns  it  rigidly  round  the  common  representative  edge  as  axis, 
leaving  its  inclination  6  intact.  Notice  also,  that  all  the  normals 
nK  being  coplanar,  the  vector  sum  in  (22)  is  again  a  vector  in  their 
plane.  Such  a  vector  can  always  be  represented  as  the  sum  of 
two  vectors  along  two  fixed  unit  vectors  in  (or  parallel  to)  that 
plane.  Let  these  be  a,  b.  Then  (22)  will  assume  the  form 

r'=r+,4a  +  £b,  (23) 

where  A,  B  are  some  scalar  functions  of  the  orientation  of  r,  and 
of  the  characteristics  of  the  prism  system,  i.e.  of  the  pK  and  the 
relative  orientation  of  the  nK  (the  latter  already  including  all  the 
angles  €K).  These  scalar  functions  A,  B  can  easily  be  determined 
for  the  case  of  any  concretely  prescribed  prism  system.  The 
resultant  deviation  (angle)  ^  =  arccos  (rr')  follows  at  once  by 
multiplying  (23)  scalarly  by  r,  which  gives 

cos  j/'  =  i  +  A  cos  (r,  a)  +  B  cos  (r,  b).  (24) 

As  a  useful  exercise  the  reader  may  treat  fully  the  case  of  two 
glass  prisms  separated  by  an  air-gap,  either  as  an  "  ideal  "  or  as 
a  "  defective  "  system. 

In  the  latter  case  it  will  be  of  particular  interest  to  fully  develop 
the  formula  for  the  effect  of  a  small  deviation  from  edge  parallelism. 

7.  System  of  Plane  Mirrors. 

Let  again,  as  in  Section  6,  each  interface  <rK  be  plane,  but  let 
each  of  them  be  contemplated  as  a  reflector.  If  the  position  of 
rays  is  disregarded,  such  a  system,  which  may  shortly  be  called 
a  multiple  reflector,  will  be  entirely  determined  by  the  assemblage, 
and  the  order  of  succession,  of  the  unit  normals  nlf  n2,  n3,  etc., 
drawn,  let  us  say,  away  from  the  reflecting  faces. 

As  was  already  hinted  at  on  p.  3,  the  resultant  formula  (IV) 
assumes  for  any  such  multiple  reflector  a  particularly  simple 
form. 

In  fact,  returning  to  the  reflection  formula  (I),  written  down 
at  the  very  beginning,  and  using  now  a  single  dash  for  every 
reflected  ray,  we  have,  for  the  first  reflected  ray, 

r^^i! -2^(^1,5,  (25) 

where  il  is  the  (original)  incident  ray.  We  already  know  that 
this  is  a  particular  case  of  the  refraction  formula  (II),  page  2, 
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viz.  with  /xr  replaced  by  —  /*.  But  there  is  this  important  difference, 
that  while  (II)  contained  in  its  right-hand  member,  through  g, 
the  "  unknown  "  s',  formula  (25)  contains  in  its  right-hand  member 
ij  only  (besides  the  given  nx)  and  no  trace  of  r^.  It  is  true  that, 
as  was  promised  on  p.  2,  the  vector  s'  can  easily  be  eliminated 
from  the  right-hand  member  ;  but  in  doing  'so  we  obtain  for  the 
refracted  ray  a  formula  which  in  its  structure  differs  essentially 
from  the  reflection  formula  (25).  In  fact,  squaring  (II),  and 
remembering  that  s/2=//2,  S2=//.2,  we  have 


whence  g,  in  terms  of  sn  =/*  cos  i,  and  by  substitution  in  (II), 

r  '  =  pi  -  /on  (nr)  +  n  V^2  +  p2  (nr)2,  (26) 

where  />=/*//*'  and  <r2  =  i  -p2.  Here  r  only  appears  in  the  right- 
hand  member  ;  but  the  operation  to  be  performed  upon  r  in  order 
to  obtain  the  refracted  ray  r'  is  *  manifestly  a  w0w-distributive  one. 
That  is  to  say,  if  (26)  is  written  for  the  moment 

r'  =  wr, 

w  standing  for  an  operator,  then  co(A+B)  is  not  equal  wA+wB. 
It  is  for  this  reason  that  we  have  said  (Section  i)  that  it  is  not 
worth  while  to  eliminate  s'  from  g.  Nor  would  it  give  us  any 
advantage  to  study  the  complicated  and  non-distributive  operator 
o>,  as  defined  by  (26),  —  which  might  be  called  the  "  refractor." 

On  the  other  hand,  a  glance  at  the  reflection  formula  (25)  will 
suffice  to  show  that  the  operator  converting  the  incident  into  the 
reflected  ray  is  distributive.  This  is  the  reason  of  the  simplicity  of 
multiple  (plane)  reflection,  corresponding  to  a  repeated  application 
of  such  operators. 

In  fact,  n^A+B)  is  njA+iijB,  which  proves  the  distributivity 
of  the  operator  fij,  say,  involved  in  (25).  That  is  to  say,  if  we 
write  that  equation  rj^O^,  we  have,  for  any  A,  B, 


This  is  a  capital  property  of  the  "  reflector  "  ft,  as  contrasted  with 
the  non-distributive  property  of  the  "refractor"  w. 

In  order  to  see  the  structure  of  the  "  reflector,"  belonging  to  the 
mirror  whose  normal  is  n1;  write  (25),  using  the  dot  (instead  of 
brackets)  as  separator, 

r1/=r1-2n1.n1r1. 

*  Quite  apart  from  being  more  complicate  than  the  reflecting  operation. 

c 
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This  is  but  a  slight  change  in  notation.  Next,  instead  of  writing 
twice  the  same  operand  rlf  put 

to  mean,  by  convention,  precisely  the  same  thing  as  before.  Thus, 
the  operator  fllf  converting  any  incident  ray  into  the  corresponding 
reflected  ray,  is  seen  to  be 

H!  being  the  mirror's  unit  normal.  Such  an  operator  is  called  a 
dyadic  (nx .  nx  being  called  a  dyad,  since  it  contains  two  vectors 
separated  by  a  dot).  The  reader  need  not,  however,  be  deterred 
by  new  names  or  new  concepts.*  Since  he  knows  how  12X  was 
introduced,  there  can  be  no  mystery  about  the  meaning  or  the 
properties  of  this  operator.  In  fact,  he  can  read  the  definition  (27) 
in  plain  English  thus  :  to  operate  with  12j  upon  any  given  vector  A 
means  to  cut  off  twice  its  projection  upon  nx  along  nx  and  to  subtract 
the  vector  thus  obtained  from  I  .  A,  that  is,  from  A  itself.  This 
is  a  very  simple  operation  indeed,  more  simple  than  most  of  those 
the  reader  daily  performs,  whether  intellectually  or  physically. 
There  will  thus  be  no  danger  of  a  misunderstanding  in  using  this 
operator,  which  is  but  a  very  particular  example  of  the  powerful 
dyadic. 

Thus  our  first  reflected  ray  will  be 

r^IVi-  (28) 

Let  now  the  reflected  ray  r/  impinge  upon  the  second  mirror  of  the 
system.  Then,  !22  =  i-2n2.n2  being  the  "reflector"  of  this 
mirror,  the  second  reflected  ray  will  be,  exactly  as  before, 

and  substituting  r/  from  (28), 

r2'=QJtiiTi,  (280) 

which  means  :  operate  with  i\  upon  rlf  then  with  fi2  upon  the 
result  of  the  first  operation.  This  is  in  essence  the  same  thing  as 
to  multiply  a  given  number  x  first  by  3  and  then  the  result  by  5. 
Both  kinds  are  intelligible  and  reasonably  simple  operations.  The 
only  important  difference  between  them  is  that  while  the  order 
of  the  operations  in  the  latter  case  is  a  matter  of  indifference,  it 

*  I  mean,  new  to  the  reader.     The  name  and  the  concepts,  due  to  Gibbs,  are 
old  enough. 
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is  not  so  in  the  former  case,  12^  is  (in  general)  not  the  same 
operation  as  £l-£l%.  So  also  is  the  actual  reflection  from  mirror  i. 
followed  by  a  reflection  from  mirror  2.  not  the  same  thing  as  a  double 
reflection  in  the  reversed  order  (unless  the  two  mirrors  are  in  a 
particular  relation  to  one  another). 

By  (2$a)  it  can  now  be  said  that  the  "  reflector  "  of  a  double 
mirror  (nl7  n2)  is 


Similarly  for  three,  four  and  any  number  of  plane  mirrors. 
The  "  reflector  "  of  a  multiple  mirror,  consisting  of  N  plane  mirrors, 
acting  in  the  order  I,  2,  ...  K,  is 


or,  more  shortly, 

12  =  1112,,,     where  12K  =  I  -  2nK  .  nK,  (29) 

the  "  product  "  extending  oVer  all  the  component  mirrors  in  their 
prescribed  order.  If  any  ray  r  impinges  upon  the  multiple  mirror, 
then  the  final  ray,  reflected  in  the  said  order,  is 

r'=12r.  (30) 

The  problem  is  thus  reduced  to  building  up  the  resultant  operator 
12  from  the  component  ones,  12K.  Now,  this  is  done  almost  as 
simply  as  the  algebraic  multiplication  of  a  succession  of  poly- 
nomials. For  to  any  succession  of  operators  12X  belongs  the 
associative  property,  i.e. 


and  similarly  for  any  number  of  component  reflectors.*  The 
said  simplicity  of  handling  reflectors  12K  and  of  building  them 
up  into  the  resultant  reflector  12  for  the  multiple  mirror  is  based 
upon  these  two  properties,  the  distributivity  and  the  associativity 
of  these  operators.  In  each  case  the  process  will  consist  in  simply 
"  multiplying  "  out  the  dyads  n  .  n  contained  in  the  component 
reflectors.  The  only  precaution  needed  is  to  keep  the  order  intact 
and  to  remember  that  juxtaposed  vectors  not  separated  by  a  dot 
are  fused  into  ordinary  scalar  products.  Thus,  for  example, 

nt  .  n^g  .  n2=n1(n1n2)  .  n2  =  «12n1  .  n2, 

where  a12=n1n2  =  cos  (n1?n2)  is  an  ordinary  scalar  number.    Similarly 

*  Henceforth  "reflector"  will  be  used  for  the  operator  and  "mirror"  for  the 
reflecting  surface  to  which  it  belongs. 
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for  the  product  of  any  number  of  dyads.     Such  a  "  product  "  will 
always  be  again  a  dyad,  including  an  ordinary  scalar  factor. 

In  all  the  following  developments  we  shall  employ  the  short 
notation 

ntnK  =  cos(nt,  nK)=a«,  (31) 

so  that,  obviously, 

aiK  =  aKL. 

Before  passing  on  to  examples,  notice  that  equation  (25)  or  (28), 
when  squared,  gives 


Thus,  the  incident  it  being  a  unit  vector,  so  is  also  the  reflected  r/. 
The  same  property  holds,  of  course,  for  any  multiple  mirror, 
i.e.  if  r'=12r,  as  in  (30),  then  r'2  =  rz  =  i.  Notice  also  that  the 
iteration  of  any  12K  is  I  or  an  idem  factor  ;  in  fact, 

^A  =  i  -4n*  .  n<  +4(nK  .  nK)2=  i  -4&K  .  nK  +nK  .  nx 
(since  nK2=  i),  whence  12A  or  fi«2=  I,  which  is  an  obvious  property. 
If  \//  be  the  angle  contained  between  the  incident  and  the  finally 
reflected  ray,  we  have 

cos^=rflr,  (32) 

for  any  multiple  mirror.  The  expression  on  the  right  hand  is  the 
scalar  product  of  r  into  the  vector  12r,  or,  which  turns  out  to  be 
precisely  the  same  thing,  of  the  vector  ri2  into  r.  12  being  a  dyadic, 
composed  of  dyads  a  .  b,  can  manifestly  be  used  either  as  a  post- 
factor  or  a  prefactor.  For  any  of  the  component  dyads  can  be  used 

so.     Thus  • 

a  .  br  means  a(br) 

and  ra  .  b  means  b(ar), 

both  being  perfectly  definite  results  (vectors),  in  general  differing 
from  one  another,  of  course.  But  r(12r)  is  the  same  thing  as  (r!2)r, 
so  that  in  (32)  and  all  similar  equations  no  brackets  are  needed. 

The  reflector  12  belonging  to  any  Af-uple  mirror  will  be  obtained 
by  multiplying  out  (with  the  said  precautions)  all  the  component 
reflectors  12*.  Thus,  to  obtain  12  for  a  quadruple  mirror,  say,  we 
must  necessarily  have  first  the  operators  for  a  triple,  and,  before 
that,  for  a  double  mirror.  Similarly  for  any  N.  It  is  therefore 
but  natural  to  begin  with  the  double  mirror,  and,  before  passing 
on  to  the  triple  and  more  complicated  mirrors,  to  stop  a  while  at 
each  of  the  intermediate  stages,  discussing  and  interpreting  the 
results.  It  will  be  understood  that  by  an  "  Af-uple  "  not  neces- 
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sarily  a  system  of  N  different  simple  mirrors  is  meant,  but  only  that 
N  successive  reflections  (in  a  prescribed  order)  are  contemplated. 
An  incident  ray  can  be  reflected  from  two  mirrors,  for  instance, 
three  or  very  many  times,  as  is  the  case,  for  example,  in  the  plan- 
parallel  piece  of  a  Lummer  plate. 

Double  Mirror.  The  cosine  of  the  angle  between  the  two  com- 
ponent mirrors  being  a]2  =  n1n2,  as  in  (31),  we  have,  by  (29),  for 
the  reflector  of  our  double  mirror, 

fi=fi2121  =  [i  -2n2  .n2][i  -21^  .  nj, 

i.e.,  with  due  regard  to  the  prescriptions  of  the  "  multiplication  " 
rule, 

ft  =  i  -  2nx .  nx  -  2n2  .  n2  +  4a12n2 .  nlf  (33) 

or,  introducing  the  vector  m  =n:  -  2#12n2, 

ft  =  i  -2m  .  nj  -2n2  .  n2.  (33#) 

This,  of  course,  is  valid  for  any  operand,  that  is  to  say,  for  any 
incident  ray  r,  giving  r'  =^1  for  the  reflected  ray.  Thus 

r/=r-2r1m-2r2n2, 

where  r1=rn1,  r2=rn2  are  the  given  projections  of  the  incident 
ray  upon  the  mirror  normals.  Since  the  edge  e,  along  which  the 
two  mirrors  intersect,  is  perpendicular  to  nx,  n2,  and  therefore  also 
to  m,  we  have  again,  as  in  the  case  of  a  prism, 


In  fact,  the  double  mirror  is  but  a  sub-case  of  the  prism  (for  /// 
replaced  by  — /*).  The  angle  9  which  the  ray  makes  with  the  edge 
remains  intact  after  the  double  reflection.  This  is  a  general 
property  of  the  double  mirror  independent  of  the  value  of  al2. 

From  (33)  we  see  that  figi^  differs  from  fi^g  since  the  last 
dyad  n2 .  n^  is  not  symmetrical  (n2,  na  are  two  different  directions). 
Thus,  if  a  beam  of  parallel  rays  r  is  broad  enough  to  impinge  upon 
both  mirrors,  it  is  split  by  the  double  mirror  into  two  beams,  r'  as 
before  and  T"  =  ^11^2T,  so  that 

r'-r"  =  4<212[n2  .  r^-ii! .  n2]r.  (34) 

In  particular,  if  the  double  mirror  is  orthogonal,  that  is  to  say, 
if  a12  =  n1n2  =  o,  we  have  r'=r" 

i.e.  the  two  reflected  beams  coincide.     In  this  case  the  operator  (33) 
becomes  a  symmetrical  (or  self -conjugate)  dyadic,  viz. 
=12^2  =  i  -  2nx .  nt  -  2n2 .  n2. 
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Again,  since  n^n^  we  have 

n1.n1+n2.n24-e.e  =  i, 

for  if  i,  j,  k,  be  any  normal  system  of  unit  vectors  (as  is  HJ,  n2,  e) 
we  have  manifestly  [i  .  i+j  .  j  +k  .  k]V=V,  for  any  vector  V. 
Therefore,  for  an  orthogonal  double  mirror,  whose  unit  edge  is  e, 

12=  -i  +2e.  e.  (35) 

That  is  to  say,  the  reflection  from  such  a  mirror  is,  independently 
of  the  order  (12  or  21),  wholly  equivalent  to  the  reflection  from  a 
simple  mirror  whose  normal  is  e,  followed  by  a  simple  reversal  (  -  i). 
This  holds  for  any  incident  ray  r.  If  the  incident  ray  happens  to 
be  normal  to  the  edge,  i.e.  if  re=o,  we  have 

r'=flr=-r, 

that  is  to  say,  the  ray  is  sent  back  parallel  to  itself.  This  simple 
property  is  also  familiar  from  ordinary  geometrical  constructions. 
But  none  of  the  cases  of  a  skew  r  (oblique  to  e)  will  be  found  so 
easy  to  treat  by  the  usual  methods,  based  on  a  good  drawing  and 
spherical  trigonometry,  as  by  the  proposed  vectorial  or  dyadical 
method. 

Returning  once  more  to  the  general  formula  (33),  we  have  for  the 
angle  ^  between  the  incident  and  the  finally  reflected  ray,  by  (32), 

i(i  -  cos  $  =  r^  +  r22  -  2al2r1r2  ; 
on  the  other  hand,  since  r=r1n1+r2n2  +  (re)e  is  a  unit  vector, 

f]2  +  r22  +  2a^2r1r2  =  I  -  (re)2  =  sin2  (9, 
so  that  the  last  equation  can  be  written 
i  -  cos      =  sin2<9  = 


giving  the  angle  ^  =  r,  r'  for  any  double  mirror  and  for  any  given 
incident  ray  r  =r1n1  +r,p.2. 

From  the  general  formula  (36)  we  see  without  difficulty  that  there 
is  no  such  double  mirror  which  would  send  back  (parallel  to  its 
own  path)  every  incident  ray  ;  in  short,  that  a  double  mirror  cannot 
be  a  "  central  "  mirror,  as  are  those  used  in  modern  signalling. 
In  fact,  cos^=  -  I  would  imply  4a12r1r2  +  (re)2  =  O,  and  this  cannot 
be  satisfied  for  all  directions  of  r. 

Further  discussions  of  the  double  mirror  will  be  left  to  the  reader 
as  an  exercise.  In  this  case  the  adaptation  of  the  vector  formulae 
to  numerical  purposes  is  almost  immediate,  and  calls,  therefore,  for 
no  further  explanations. 
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Triple  Mirror.  The  resultant  reflector  in  this  case  is  fi  ^ 
the  product  of  (23=i  -2n3  .  n3  into  the  operator  of  the  preceding 
case  already  developed  in  (33).  Thus,  for  any  triple  mirror,  the 
order  of  reflections  being  I,  2,  3,  the  resultant  operator  12  =^123  1S 

fl123  =  i  -  2[nj  .  nx  +n2  .  n2  +n3  .  n3] 

.  nj  -  8012tf23n3  .  nr    (37)* 


Here  again,  the  second  term  being  symmetrical,  but  not  the  third, 
nor  the  last,  12132r  etc.,  will  differ  from  ft123.  In  other  words,  an 
incident  beam  of  parallel  rays  r  will  give  rise  to  six  reflected  beams 
represented  by 


which  will,  generally  speaking,  deviate  from  one  another.  These 
six  reflected  beams  become  parallel  to  one  another,  i.e.Q  becomes 
independent  of  the  order  of  the  three  successive  reflections  when 
and  only  when 

*12  =  ^23  =  031  =  °> 

abolishing  the  third  and  the  fourth  terms  in  (37),  i.e.  when  the 
three  component  mirrors  are  mutually  perpendicular.  Moreover, 
n1}  n2,  n3  being  in  that  case  a  triad  of  normal  unit  vectors,  the 
dyadic  ^  .  n!+n2  .  n2+n3  .  n3  becomes  an  idemf  actor  or  I,  and 
the  operator  is  reduced  to 

fi«-l;     .'.  r'=-r.  (37*) 

That  is  to  say,  every  incident  ray,  independently  of  the  order  of  its 
three  reflections,  is  sent  back  parallel  to  itself.  The  orthogonal 
triple  mirror  is  a  central  mirror.  Such  triple  mirrors,  in  fact,  are 
used  in  practice,  their  orthogonality  being  now  attainable  to 
a  high  degree  of  approximation, 

Returning  to  the  general  triple  mirror,  write  rlt  r2,  r3  for  the 
direction  cosines  of  the  incident  ray  or  its  projections  upon  the 
three  normals  n*.  Then  the  direction  cosines  of  the  reflected 
ray  r'  will  be  linear  homogeneous  functions  of  rlt  r2,  rs,  viz.,  by  (37), 


and  so  on. 

*  This  may,  at  first  sight,  seem  very  complicated-  Try,  however,  to  treat  the 
general  triple  mirror  by  the  ordinary  methods.  Then  (37)  will  be  found  very 
simple. 
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Exercises.  Write  down,  by  the  general  formula  (37),  the  remaining 
two  direction  cosines  f8',  ra',  thus  obtaining  three  scalar  formulae 
with  DJ,  n2,  n3  as  reference  system. 

Assuming  that  nx,  nt,  D3  are  not  coplanar,  i.e.  that  the  three 
mirrors  constitute  a  proper  pyramid  (not  a  prism),  any  incident  ray  r 
can  be  put  into  the  form 


where  plt  etc.,  are  some  scalars.     Determine  them  in  terms  of  the 
orthogonal  projections  f1=rn1,  etc.     (Notice  that  pl=r1>  etc.,  when 
and  only  when  the  mirrors  are  mutually  perpendicular.) 
Similarly  determine  the  coefficients  of 


and  replace  the  formulae  (38)  by  three  equations  giving  the  pK'  in 
terms  of  the  pK,  and  in  terms  of  the  properties  of  the  mirror,  i.e.  a23, 
etc.  Discuss  the  general  formulae  thus  obtained,  and  apply  them  to 
a  number  of  particular  cases. 

Treat  the  prismatic  arrangement  of  the  three  component  mirrors, 
when  nlf  n2,  na  are  coplanar  ;  i.e.  such  that  an^  +6n2  +  cn3  =  o,  a,  b,  c 
being  three  scalars,  and  so  on.  Consider  first  an  incident  ray  j_  to 
the  edges  (all  represen  table  by  one  e),  and  then  any  oblique  or  skew 
ray.  In  the  former  case  imagine  the  mirror  pierced  at  a  place,  just 
to  let  through  the  incident  ray  ;  or  else  let  the  mirrors  consist  of 
non-silvered  (plane  parallel)  glass  plates. 

The  angle  ^  between  the  incident  and  the  finally  reflected  ray 
is  given  by  cos  i/^rfir,  i.e.,  according  to  (37),  by 

J(i  -  cos  ^)  =  (rf  +  r£  +  r32)  -  2  [alzrtr2  +  etc.]  +  4^12^23^1-  (39) 
Each  of  the  two  bracketed  expressions  is  symmetrical  with  respect 
to  the  suffixes  I,  2,  3,  but  the  last  term  is  not.  Thus,  if  we  replace 
the  order  of  reflections  123,  by  132,  say,  and  call  ^'  the  angle 
between  the  corresponding  finally  reflected  ray  and  r,  we  have  the 
elegant  formula 

cosf-  cos  r/-  =  8r1a23(a12r3  -  al3r2)  .  (40) 

Similarly  for  any  other  permutation.  Notice  that  321  gives  the 
same  reflected  ray  as  123. 

If  the  three  mirrors  constitute  a  regular  pyramid,  i.e.  if 

a12  =  a23=asl  =  COS  WJ  SaY> 

then  formula  (39)  becomes 

!(  I  -  cos  \//)  =  (r-f  +  r22  +  rs2)  -  2  cos  w(rxr2  +  r2rs  +  r3rt) 

.  (390) 
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More  especially,  if  the  incident  ray  r  is  equally  inclined  to  the  three 
reflecting  planes  of  the  regular  pyramid,  we  have 


• 

and  therefore,  remembering  that  I  -  cos  ^  =  2  sin2  —  , 


=cot-/v/i  - 


sin-=cot-/vi  -2  cosw-f--cos2w. 


In  this  case  the  angle  \ff  is,  of  course,  independent  of  the  order  of 
reflections.  The  reflected  beams  (corresponding  to  a  beam  of 
parallel  incident  rays)  are  not  parallel  to  one  another,  but  are 
equally  inclined  to  and  symmetrically  disposed  around  the  direction 
of  the  incident  beam.  These  reflected  beams  will  coincide  in 
direction  when  and  only  when  w  =  90°,  i.e.  when  the  mirror  becomes 
an  orthogonal  and,  therefore,  a  central  mirror.  In  fact,  for 

w=-,  formula  (396)  gives  sin-  =  i,  i.e.  ^  =  180°,  which  means  a 

reversal  of  the  incident  ray.  This,  as  we  already  know,  is  the 
property  of  an  orthogonal  mirror  with  respect  to  any  incident  ray. 
Quadruple  and  more  Complicated  Mirrors.  Multiplying  (37)  by 
I  -  2n4  .  n4,  where  n4  is  the  normal  of  the  fourth  mirror,  the  reader 
will  obtain  the  operator  ^=^^34,  and  similarly  for  quintuple 
reflections,  and  so  on.  As  was  already  remarked,  some  of  the 
normals  n*  may  reappear,  that  is,  we  may  have,  for  instance,  a 
quintuple  reflection  from  but  two  mirrors  ;  if  such  be  the  case, 
it  is  enough  to  put  n3=n5=n1,  and  n4=n2. 

As  an  example  of  this  kind  the  reader  may  treat  the  rotation  and 
progress  of  a  skew  ray  incident  upon  an  interior  face  of  an  equilateral 
prism  built  up  of  three  plane  mirrors  ;  nlf  n2,  na  coplanar,  n4=xii, 
ns  =na,  and  so  on. 

Reversal  of  the  Order  of  Reflections.  For  any  multiple  mirror,  let 
r  be  the  incident  ray,  r'  the  finally  reflected  ray  when  the  order  of 
reflections  is  123  ...  N  -  I,  N,  and  r"  the  finally  reflected  ray  when 
the  order  of  reflections  is  reversed,  i.e.  N,  N-  1,  ...  321.  Then, 
ft  being  fiA_i  ...  fl3QaQlf  as  on  page  25,  we  have 

r'=ftr     and    r"=rfi, 
whence 

rr'=r!2r,     r"r  =rr"=rfir, 

and  therefore,  for  any  multiple  mirror, 

r'r=r"r,  (41) 
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while  r'r"=rfi2r.  That  is  to  say,  the  reflected  rays  r',  r",  although 
not  parallel  to  one  another,  are  always  equally  inclined  to  the  incident 
beam  r.  The  equality  ^123  =  ^321  mentioned  on  page  30  is  but 
a  special  instance  of  this  general  property. 

8.  Dyadic  representing  the  most  General  Reflector. 

Consider  the  successive  reflections  of  a  ray  r  from  any  number 
of  mirrors,  whose  normals  are  nlf  n2,  etc.,  and  let  ft,  with  the 
order  of  reflections  subintended,  be  the  corresponding  resultant 
operator  or  the  "  reflector "  belonging  to  the  whole,  ordered, 
system.  (This  comprises  also  the  case  of  any  number  of  curved 
mirrors  ;  the  only  difference  being  that  when  they  are  plane,  all 
the  normals  are  immediately  given,  and  when  curved,  the  successive 
normals  are  to  be  found  by  the  application  of  the  "  transfer  for- 
mula," developed  in  a  previous  Section.)  Now  ft  being  the  product 
of  dyadics  of  the  form 

I  -  2n  .  n, 

is  itself  a  dyadic.  And  since  our  space  has  but  three  dimensions, 
all  the  consequents  *  of  any  reflector,  no  matter  how  complicated 
the  system  to  which  it  belongs,  can  be  expressed  as  linear  homo- 
geneous functions  of  three  non-coplanar  vectors,  which  can  always 
be  made  unit  vectors,  let  us  say,  f,  g,  h.  The  antecedents  belonging 
to  these  consequents  will  then,  for  any  given  multiple  mirror,  be 
certain  three  vectors,  let  us  say,  a,  b,  c.  Thus,  the  most  general 
reflector  will  have  the  form  of  a  trinomial  dyadic, 

ft=a  .  f  -fb  .  g+c  .  h.  (42) 

The  plain  meaning  of  this  is  that  any  incident  ray  r  will  be  reflected, 
in  the  prescribed  order,  as 

r'  =  a(fr)  +b(gr)  +c(hr).  (420) 

The  consequents  f,  g,  h  need  not  necessarily  constitute  a  normal 
system  ;  it  may  sometimes  be  convenient  to  choose  an  oblique 
system  of  vectors.  The  choice  will  depend  upon  the  properties  of 
the  system  of  mirrors. 

The  above  result  may  also  be  stated  by  saying  that  the  most 
general  ft  is  a  linear  vector  operator,  or  r'  a  linear  homogeneous 
vector  function  of  r,  as  was  already  mentioned. 

*  According  to  Gibbs'  terminology,  the  first  vector,  a,  in  a  dyad  a .  b  is  called  the 
antecedent,  and  the  second  vector,  b,  the  consequent. 
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Splitting  into  Cartesians,  for  instance,  and  counting  up  the 
coefficients  expressing  the  components  of  r'  in  terms  of  those  of  r, 
the  reader  will  see  at  once  that  the  most  general  linear  vector 
operator  implies  nine  independent  scalar  or  numerical  data.  Let 
us  count  the  number  of  such  data,  implied  in  the  most  general 
reflector  ft,  upon  the  formula  (42)  itself. 

The  consequents,  f,  g,  h,  which  can  always  be  taken  as  unit 
vectors,  and,  moreover,  always  as  a  normal  system  of  unit  vectors,* 
imply  2  +  1+0=3  scalar  data  (for  the  direction  f  being  fixed,  g 
requires  but  one  datum,  and  h  is  codetermined,  apart  from  the 
irrelevant  sign  or  sense).  Now  it  is  always  (i.e.  for  any  trinomial 
dyadic)  possible  to  give  to  the  normal  system  f,  g,  h  such  an  orienta- 
tion as  to  convert  the  antecedents  also  into  some  normal  system  of 
vectors,!  but,  of  course,  not  unit  vectors,  after  f,  etc.,  have  been 
made  so.  Thus  the  orientation  of  the  system  of  antecedents 
implies  three  more  data,  and  the  size  or  absolute  value  of  each 
antecedent  one  datum.  Thus  the  total  number  of  scalar  data 
implied  in  the  most  general  trinomial  dyadic  is  nine. 

And  since  our  reflector  ft  is  such  that  it  leaves  the  size  of  any 
operand  unaltered,  i.e. 

r'2=r* 

for  any  r,  the  number  of  independent  scalar  coefficients  contained 
in  the  most  general  reflector  ft  is  reduced  to  eight.  These  coefficients 
depend  upon  the  properties  of  the  system  of  mirrors  and,  of  course, 
of  the  prescribed  order  of  successive  reflections.  No  matter, 
however,  what  their  arrangement  and  how  many  of  simple  mirrors 
are  built  up  into  the  system,  all  of  its  properties  are  crowded  into 
eight  coefficients, — corresponding  to  each  prescribed  order.  On 
permutating  the  order,  the  operator,  or  its  coefficients,  will,  in 
general,  be  changed  ;  in  some  special  cases  they  may  remain 
intact.  The  reversal  of  order  has,  however,  the  perfectly  general 
property  expressed  by  (41). 

A  dyadic,  or  linear  vector  operator,  such  that  ftr  =  rft,  for  any  r, 
is  called  a  self -conjugate  or  symmetrical  one.  The  simple  reflector 

*  Although  this  may  sometimes  be  not  the  most  convenient  choice.  This,  how- 
ever, does  not  influence  the  result  of  our  counting. 

t  A  simple  proof  is  given  in  J.  W.  Gibbs'  lectures  on  Vector  Analysis,  edited  by 
E.  B.  Wilson,  New  York,  1907,  p.  304.  Our  dot,  a  separator,  is  not  to  be  con- 
founded with  Gibbs5  dot,  which  is  the  symbol  of  scalar  multiplication.  The 
notation  adopted  in  the  present  volume  is  essentially  that  of  Oliver  Heaviside. 
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I  —  2n  .  n  is,  obviously,  such  an  operator.  But  that  corresponding 
to  the  general  double  mirror  is  no  longer  self  -conjugate.  Other 
details  of  the  discussion  of  the  multiple  reflection  formulae  may 
be  left  to  the  reader.  As  to  further  examples  and  exercises, 
besides  those  proposed,  the  reader  will  doubtless  invent  them  himself 
without  difficulty. 

9.  Differential  Properties  of  Refracted  or  Reflected  Pencils. 

Although  the  chief  object  of  this  book,  viz.  the  tracing  of  a  single 
ray,  is  in  its  essence  exhausted  by  what  has  been  said  in  the  pre- 
ceding sections,  it  may  be  well  to  add  one  more  section  in  order 
to  give  the  reader  at  least  some  cursory  hints  about  the  treatment 
of  differential  properties  of  beams  or  pencils  of  rays  sent  through  an 
optical  system  :  properties,  that  is  to  say,  which  concern  the 
passage  from  one  ray  of  a  pencil  to  its  neighbours,  or  the  passage 
from  a  given  wave-length,  and  the  corresponding  index  //,  to  an 
infinitesimally  different  one  (astigmatism,  chromatic  aberration). 

Problems  of  this  kind  do  not,  as  a  matter  of  fact,  require  anything 
essentially  different  from  what  has  already  been  given  in  the 
earlier  part  of  the  book.  They  are  solved,  in  each  case,  by  simply 
differentiating  both  sides  of  the  vector  formulae  concerning  a 
single  ray,  obtainable  by  the  above  method.  Now,  the  differentia- 
tion of  vector  expressions  does  not  differ  from  that  of  scalar  ones, 
provided  —  of  course  —  that  we  differentiate  with  respect  to  a 
scalar  magnitude  as  the  independent  variable  or  parameter.  As 
a  precaution  hardly  anything  more  is  needed  than  the  preservation 
of  order  of  the  (vector)  factors  of  a  vector  product  and,  in  case 
of  dyadics,  the  discrimination  between  antecedent  and  consequent. 

Thus,  if  f  means  the  "  differential  coefficient  "  or  the  flux  of  r 

c)r 
"with  respect  to  any  scalar  variable  x,  i.e.  if  r  stands  for  ^-,   the 

flux  of  a  sum  a+b  of  two  vectors  is  a+b,  and  that  of  the  scalar 
product  ab  is  ab  +ab.     Similarly,  for  the  vector  product, 


and  so  on.     Or,  using  "  differentials," 

d(&  +  b)  =  da,  +  db,     rf(ab)  =  a  d\>  +  b  rfa, 
where  da,,  db  are  infinitesimal  vectors.     Similarly 
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the  second  term  being  negatived  to  compensate  for  the  reversal  of 
the  order  of  its  two  factors.  Further  explanations,  if  at  all  needed, 
will  be  found  in  any  treatise  on  Vector  Analysis.* 

Such  being  the  case,  it  will  be  enough  briefly  to  consider  the 
passage  of  a  pencil  of  rays  r  through  a  single  refracting  surface 
(reflection  being,  formally,  but  a  sub-case).  Thus,  returning  to 
the  fundamental  refraction  formula  (II),  Section  I,  we  have  for 
any  ray  of  an  incident  pencil,  let  us  say  for  the  "  central  "  ray  of 

the  pencil, 

s'  =  s+£n,  (II) 

where  n  is  the  normal  of  the  refracting  surface,  of  any  shape  (free 
of  singularities),  at  the  point  struck  by  the  central  incident  ray  s. 
Let  the  refractive  indices  /*,  /*'  be  fixed,  and  let  us  pass  from  the 
central  ray  s  (to  which  corresponds  the  central  ray  s'  of  the  refracted 
pencil)  to  an  infinitesimally  neighbouring  ray  s  +  ds.  If  all  the 
rays  of  the  pencil  are  parallel,  then,  of  course,  ds=O,  s  representing 
the  whole  pencil.  But  in  general  we  shall  have  ds^o.  At  any 
rate,  since  s2=/>t2,  the  relation  sds=o,  or  ds±s,  will  hold.  The 
corresponding  refracted  ray  will  be  s'-t-ds',  and  ds'_Ls'  as  before. 
Now,  passing  from  one  ray  to  another  means,  in  general,  varying  the 
normal  n,  and  therefore  also  g.  Thus,  subtracting  the  original  or 
central  s'  from  the  new  vector  s'-fds',  we  have  for  the  required 
change  of  the  refracted  ray, 

ds'  =  ds+gdn+n  dg,  (43) 

where,  by  the  definition  of  g,  p.  2, 

dg=f*  sin  i .  di  —  p  sin  i'.  di't 
and  by  the  sine-law,  //  cos  i'.  di'  =  /x,  cos  i .  di,  so  that  ultimately 

dg=f(i)di, 

f  being  a  known  function  of  i  alone.  Notice  that,  n  being  a  unit 
vector,  we  have  always  ndn=o.  In  the  case  of  reflection  it  is 
enough  to  replace  //  by  — /x. 

In  essence,  formula  (43)  already  contains  the  answer  to  all 
thinkable  questions  concerning  the  properties  of  the  infinitesimal 
refracted  pencil  in  terms  of  those  of  the  incident  one.  It  is  enough 
to  read  it  intelligently  and  to  develop  it  appropriately  in  each 
particular  case. 

*  Cf.,  for  instance,  pp.  22-48  of  my  Vectorial  Mechanics.  Only  a  small  part  of 
what  is  there  given  will  be  required  for  the  purposes  of  geometrical  optics. 
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Let  u,  v  be  any  curvilinear  coordinates  denning  the  position  of 
a  point  at  the  refracting  surface  <r.  Then  "  incident  pencil  given  " 
means  'ds  given  as  a  function  of  u,  v  throughout  an  infinitesimal 
region  round  u0,  y0,  the  point  struck  by  the  central  ray.  Now,  the 
surface  o-  being  given  (round  that  place  at  least),  say,  through  the 
expression  of  its  squared  line-element, 

dl2  =  Eduz  +  2F  dudv  +  G  dv2,  (44) 

dn  is  also  known  as  a  function  of  the  two  coordinates.  The  same 
is  true  of  dg,  which  differs  from  di  by  a  known  factor  f(i). 

Thus  (43)  finds  ds'  as  function  of  u,  v  round  w0,  v0 ;  that  is,  it 
gives  all  the  properties  of  the  refracted  pencil.  For  s',  n  or  s', 
uo,  vo  g*ve  tne  direction  and  the  position  of  the  central  pencil,  and 

s'  +  ds'  with  u0  +  du,     v0  -f  dv 

localize  and  determine  the  direction  of  any  ray  of  the  refracted 
pencil.  Thus  every  question  concerning  the  latter  pencil  is  reduced 
to  purely  geometrical  constructions. 

Employing  u,  v  (for  instance,  <f>,  &  for  a  sphere)  as  our  independent 
variables,  we  have,  of  course, 


(du         'dv 
and  in  quite  the  same  way 


,  ,  , 

dn=^-du+^-  dv. 
ou         ov 

If  the  incident  pencil  consists,  for  instance,  of  parallel  rays,  so 
that  ds=o,  then 

dg,  (45) 


whence  nds'  =  dg,  and  Vn  ds'  =#Vn  dn,  two  convenient  formulae 
replacing  (45).  The  chief  question  in  any  such  investigation  is 
whether  the  neighbour  ray  passing  through  P(uQ  +  du,  v0  +  dv] 
meets  at  all,  at  finite  or  infinite  distance,  the  central  refracted  ray 
passing  through  P0(w0,  v0).  In  order  to  put  this  question  in 
vectorial  language,  introduce  the  infinitesimal  vector  dl,  drawn  on 
the  surface  from  P0  to  P,  its  size  dl  being  determined  by  the  quad- 
ratic form  (44)  with  E,  F,  G  as  given  functions  of  u,  v.  Then  the 
said  rays  meet  or  not  (are  coplanar  or  not)  according  as  the  scalar 
magnitude 
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vanishes   or  not.     Remembering  that  Vs's'=0,   identically,   this 
criterion-expression  is  reduced  to 


which,  by  (II)  and  (45),  is 

Q  =  [dg  Vsn  +  gVs  dn  +  #2Vn  dn]  dl.  (46) 

Given  the  refracting  surface  and  the  incident  pencil  s,  together 
with  the  central  normal  n,  this  expression  will  vanish  or  not  accord- 
ing to  the  choice  of  the  point  P,  i.e.  of  the  vectorial  line  element  dl. 
In  general  we  shall  have  Q  =o  only  for  certain  particular  directions 
of  dl  upon  the  interface.  Having  established  the  conditions 
for  Q  =  o,  it  will  be  a  comparatively  easy  task  to  localize  the  point 
of  intersection  of  the  P-ray  with  the  central  ray  of  the  refracted 
pencil,  and  thus  to  build  up  the  caustics.  If  the  incident  beam 
is  normal,  then,  whatever  the  shape  of  the  refracting  surface, 
the  most  convenient  curvilinear  coordinates  u,  v  will  be  those 
measured  along  the  lines  of  curvature  of  that  surface  drawn  through 
P0  as  origin,  in  which  case  (44)  is  reduced  to  dl2  =  E  duz  +  G  dvz. 
For  oblique  incidence,  however,  another  network  may  in  certain 
cases  be  found  more  convenient. 

As  an  exercise  the  reader  may  treat  the  case  of  normal  incidence 
upon  any  refracting  surface,  such  as  an  ellipsoid,  a  cylinder  or 
a  cone  ;  and  the  case  of  oblique  incidence  for  a  spherical  surface, 
thus  deriving  vectorially  the  well-known  formulae  for  the  sagittal 
and  the  meridional  section  of  a  pencil. 

Since  these  and  similar  problems  have  nothing  to  do  with  "  ray- 
tracing,"  the  proper  argument  of  this  little  volume,  we  cannot 
enter  upon  them  any  further.  The  last  mentioned,  and  similar, 
exercises  may,  however,  be  warmly  recommended  to  the  reader 
as  good  and  interesting  opportunities  of  training  in  the  use  of 
vector  language. 
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